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ABSTRACT 


Measurements  of  temperature  Increases  resulting  from  expo- 
sures to  argon  laser  Irradiations  were  made  In  the  eyes  of  living 
rhesus  monkeys  with  10-20vm  diameter  thermocouples.  Using  the 
appearance  of  an  ophthalmoscoplcally  visible  lesion  as  the  criterion 
for  damage,  the  temperature  associated  with  the  appearance  of  a 
threshold  lesion  5 minutes  post  exposure  was  determined.  Tempera- 
ture Increases  at  the  center  of  the  lesion  and  at  the  radial  extent 
of  the  lesion  were  determined  from  radial  scans  of  the  thermocouple 
through  the  laser  Image  and  from  ophthalmoscopic  measurements  of  the 
lesion  radius. 


Temperatures  associated  with  threshold  lesions  were 
measured  for  a number  of  exposure  durations  and  for  image  sizes 
from  100-200um  (half-power  diameter).  Thresholds  were  acquired 
for  both  macular  and  paramacular  (temporal)  exposure  sites.  The 
average  macular  threshold  temperature  increases  at  beam  center 
were  66.1°C,  29.4°C,  24.1°C  and  19.9°C  for  exposure  durations  of 


.01  second.  .1  second.  1 second,  and  10  seconds,  respectively. 
Average  paramacular  threshold  temperature  rises  at  the  center  of 
the  lesion  were  55.9°C,  39.9°C,  28.7°C  and  24.8°C  for  .01  second, 
.1  second,  1 second,  and  10  seconds  exposure  durations. 

The  empirically  determined  threshold  temperatures  were 
compared  to  a finite  differences  model  solution  of  the  heat 


conduction  equation.  The  threshold  corneal  power,  measured  retinal 
Image  profile  and  exposure  duration  were  Input  to  the  model  for 
threshold  temperature  prediction.  The  parameters  of  the  model  were 
adjusted  to  provide  agreement  between  measured  and  calculated  results 
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CHAPTER  I 


INTRODUCTION 


Intense  light  sources  such  as  Xenon  lamps  and  lasers 
are  used  In  a variety  of  applications.  Industrial  and  medical 
applications  Include  welding,  cutting,  ranging  and  photocoagula- 
tion  of  tissues.  In  any  application  of  Intense  light  sources, 
a major  consideration  Is  the  safe  use  of  the  source  In  the 
particular  application.  The  accidental  overexposure  of  personnel 
must  be  prevented.  The  characteristics  of  the  source  and  the 
mechanisms  of  Interaction  of  Intense  light  with  biologic  tissues 
are,  therefore,  very  Important. 

The  skin  and  the  eye  are  most  vulnerable  to  visible, 
ultraviolet  and  Infrared  sources.  The  ability  of  the  eye  to 
focus  visible  light  to  extremely  small  spot  sizes  makes  the  eye 
the  most  vulnerable  biologic  tissue.  The  unique  structure  of  the 
retina  and  the  fundus  of  the  eye  make  these  structures  extremely 
sensitive  to  Intense  light  energy  since  the  pigmentation  In  the 
pigment  epithelium  (P.E.)  and  choroid  absorbs  a large  amount 
of  the  light  reaching  these  layers.  Heat  conduction  from  these 
layers  damages  the  neural  retina  which  has  less  ability  for 
repair  than  skin. 

Past  studies  of  damage  to  the  fundus  of  the  eye  from 
Xenon  and  laser  radiation  Indicate  that  the  major  damage  mechanism 


!s 


Is  due  to  temperature  Increase  In  the  tissue  (1,  2,  3),  es- 
pecially for  exposure  durations  longer  than  a few  microseconds. 
For  shorter  duration  exposures,  different  physical  mechanisms 
such  as  acoustic  transients  and  steam  production  may  be  Important 
(1,  2,  3).  For  purposes  of  determination  of  safety  standards 
for  laser  radiation  and  the  extent  damage  In  clinical  uses  of  the 
Intense  light  sources,  the  determination  of  temperature  Increase 
In  the  tissue  for  a given  source  and  power  level  must  be  known. 

In  this  research,  the  primary  objective  was  to  obtain 
measurements  of  the  temperature  field  for  various  exposure  dura- 
tions with  thin  film  copper-nlckle  thermocouples  In  the  fundus  of 
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CHAPTER  II 
BACKGROUND 
Threshold  Studies 


Several  excellent  reviews  on  threshold  bum  studies 
have  appeared  In  the  literature  (1,  2,  3).  This  section  will 
explain  the  factors  In  the  study  of  threshold  bums  which  are 
pertinent  to  this  research. 

Early  investigations  of  damage  to  the  eye  as  a result 
of  exposure  to  laser  sources  were  carried  out  by  Zaret  (4)  and 
Ham  et  al.  (5).  Ham's  experiments  with  ruby  lasers  and  white 
light  sources  have  been  the  basis  for  defining  most  of  the  "safe" 
levels  currently  In  use.  The  criterion  for  damage  was  the  ap- 
pearance of  an  ophthalmoscoplcally  visible  lesion  at  times  of  5 
minutes  and  one  hour  post  exposure.  Vassllladis  et  al.  (6)  at 
Stanford  Research  Institute  found  thresholds  four  times  those 
reported  by  Ham  et  al.  (5).  Other  Investigators  have  found 
thresholds  for  various  wavelengths  and  exposure  durations  which 
vary  by  more  than  an  order  of  magnitude  from  the  threshold  values 
published  by  Ham  and  Vassllladis  (5*  6). 

Most  thresholds  have  been  characterized  by  the  power 
entering  the  pupil  or  by  the  retinal  Irradlance  In  W/cm2  or 
J/cm2,  respectively.  The  data  of  Ham  et  al.  (5)  Indicate  that 
a higher  power  density  Is  required  to  produce  a threshold  lesion 
for  a 10ym  Image  than  for  a 800ym  Image  at  long  times  of  exposure. 


For  times  of  exposure  below  one  millisecond  the  data  Indicate 
the  same  power  density  Is  required  to  produce  a lesion  for  all 
Image  sizes  for  a given  exposure  duration. 

In  Figure  1 (3),  a constant  energy  "density"  ratio 
between  thresholds  determined  by  different  Investigators  for 
Image  sizes  of  100pm  to  800pm  Is  Indicated.  The  data  Indicate 
a higher  retinal  dose  (J/cm2)  for  a smaller  Image  than  for  a 
larger  Image  for  all  times  of  exposure.  If  Injury  Is  a function 
of  retinal  temperature,  one  would  expect  a smaller  Image  to 
cause  a lower  temperature  Increase  for  the  same  value  of  retinal 
Irradlance  than  a larger  Image  because  of  conduction  effects. 

Variations  In  thresholds  determined  by  different 
Investigators  are  due  to  several  factors.  The  threshold  values 
depend  most  upon  the  spectral  and  temporal  characteristics  of  the 
source.  In  the  visible  to  near  Infrared  spectra,  most  of  the 
light  Incident  upon  the  cornea  is  transmitted  to  the  retina 
(7,  8,  9,  10).  In  the  far  Infrared  spectrum,  most  of  the  energy 
Is  absorbed  In  the  cornea.  For  example,  Vassllladls  et  al.  (6) 
report  a ratio  of  Neodymium  (Nd)  to  ruby  laser  threshold  energy 
of  approximately  five  to  one.  The  Neodymium  wavelength  (1.06pm) 
falls  In  a water  absorption  band  In  the  near  Infrared.  For  this 
reason,  more  energy  Is  absorbed  In  the  preretlnal  ocular  media 
than  for  wavelengths  In  the  visible  region. 


EXPOSURE  DURATION  IN  SECONDS 


FIGURE  1 


Retinal  Injury  threshold  energy  density 
versus  exposure  duration.  Upoer  curve 
Is  from  Vassllladls  et  al.  (6)  for  100pm 
Images.  Lower  curve  Is  from  data  of  Ham 
et  al.  (5)  for  800pm  Images. 
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Another  factor  which  affects  the  threshold  value  Is 


the  variation  In  thresholds  between  experimental  species  and 


variations  In  pigmentation  within  species.  Some  Investigators 


use  monkeys,  others  rabbits,  and  a few  studies  have  been  performed 


with  humans  (11).  Rabbit  thresholds  can  vary  by  a factor  of 


4 to  1 for  light  to  dark  pigmentation  In  rabbits  (1).  In  the 


monkey,  the  macular  and  paramacular  regions  exhibit  differences 


In  pigmentation  and  blood  supply.  The  macular  area  of  the  eye 


Is  darker  and  has  less  blood  supply  In  the  neural  layers  than 


the  paramacular  portion  of  the  fundus.  Rhesus  monkey  thresholds 


for  short  times  of  exposure  In  the  macular  and  paramacular  (nasal) 


regions  of  the  eye  differ  by  a factor  of  greater  than  2 to  1 (12). 


For  long  times  of  exposure  the  macular  and  paramacular  thresholds 
are  within  20%.  Data  from  Technology,  Inc.,  (13)  (T.I.)  Indicate 
a nearly  constant  30%  difference  between  macular  and  paramacular 


(temporal)  thresholds. 


Threshold  values  also  differ  because  of  differences  In 


lesion  grading  criteria  between  Investigators.  For  times  of 
observation  of  one  hour  or  longer,  the  threshold  for  an  ophthal- 
moscoplcally  visible  lesion  decreases  relative  to  the  5 minute 
observation.  Davis  and  Mautner  (14)  report  a value  of  80%  of 


the  threshold  determined  for  an  Inmedlately  observable  criterion 


for  an  observation  time  of  24  hours;  the  10  minute  obsqgyable 
lesion  threshold  was  90%  of  the  value  obtained  for  an  Immediately 
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observable  lesion.  At  T.I.  (13)  a time  of  5 to  7 hours  post 
exposure  has  been  found  to  be  the  minimum  observation  time 
required  to  define  a minimal  lesion.  Gibbons  (15)  finds  that 
the  one  hour  and  24  hour  thresholds  are  approximately  equal  for 
exposure  durations  shorter  than  5 seconds.  For  a 5 second 
exposure,  his  experiments  Indicate  a 5%  decrease  In  the  threshold 
corneal  power  required  to  produce  a visible  lesion.  For  120 
seconds  exposure  duration,  the  24  hour  threshold  Is  88%  below 
the  one  hour  threshold. 

In  addition  to  the  experimental  differences  mentioned 
above,  which  give  threshold  variations  of  more  than  an  order  of 
magnitude,  several  physiologic  factors  may  enter  into  the 
threshold  determination.  These  factors  include  pupil  diameter, 
body  temperature,  accommodati on  of  the  eye,  blood  pressure  and 

1 1 

level  of  anesthesia.  The  pupil  diameter  is  important  in  the 
evaluation  of  hazards  from  diffuse  reflections  or  large  beams 
for  which  all  the  energy  In  the  beam  may  not  enter  the  eye. 

The  level  of  accommodation  of  the  eye  determines  the  total  focusing 

I ■ 
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power  of  the  ocular  system  and  affects  the  minimal  image  size  for  ! 

both  "near"  and  "far"  field  hazards.  In  the  "far  field"  case 

the  eye  is  relaxed  (focused  at  Infinity)  and  not  acconmodated 

at  the  laser  as  In  the  near  field  case.  The  far  field  case  yields 

the  minimal  Image  size  for  laser  exposure  to  the  eye.  The 

theoretical  relation  between  retinal  Irradiance  and  source 
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radiance  for  an  extended  source  is  obtained  using  the  Airy 
formula  (16).  The  retinal  Irradlance  Ir  Is  given  by  the 
relationship: 
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where  de  Is  the  pupil  diameter  In  cm,  R Is  the  source  radiance 
(W/cm^  ste radian),  T Is  the  transmission  of  the  ocular  media, 
and  f Is  the  focal  length  of  the  eye.  The  formula  above  Is  for 
viewing  a diffuse  specular  reflection.  The  theoretical  lower 
limit  of  the  retinal  Image  diameter  for  the  Airy  disc  Is  given 
by: 


dr  * 


2.44X  f 


where  de  and  f are  defined  as  above  and  X Is  the  radiation  wave- 
length In  cm.  From  these  considerations,  the  dark  adapted 
(dlalated  pupil)  viewing  situation  with  the  eye  focused  at 
Infinity  presents  the  "worst"  case  for  viewing  a laser  source. 

Most  experiments  have  been  performed  under  the  worst 
case  situations  mentioned  above.  This  undoubtedly  results  In 
a much  lower  threshold  than  In  more  realistic  exposure  situations 
which  would  occur  In  cases  of  accidental  exposure.  The  subject 
may  or  may  not  be  dark  adapted  as  In  the  majority  of  laboratory 
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and  eliminates  the  blink  reflex.  Long  exposure  durations  could 


hardly,  therefore,  be  expected  to  present  a significant  hazard 
to  human  subjects  when  the  exposure  Is  near  the  threshold  value 
determined  under  the  Idealized  situations  mentioned  above. 

The  role  of  body  temperature  In  the  determination  of 
the  ophthalmoscoplcally  visible  threshold  has  been  reported  by 
Ward  and  Bruce  (17),  and  Polhamus  and  Welch  (18).  Ward  and  Bruce 
found  an  Increase  In  retinal  Irradlance  for  a threshold  lesion 
at  a lower  core  temperature  and  a decrease  In  threshold  Irradlance 


for  higher  core  temperatures.  Ward  and  Bruce  extrapolated 
their  data  to  the  rectal  temperature  axis  on  a plot  of  threshold 
Irradlance  versus  rectal  temperature  and  concluded  that  a fundus 
temperature  Increase  of  6 to  7°C  would  be  sufficient  to  cause 
damage. 

Polhamus  and  Welch  (18)  measured  the  temperature  rise 


required  to  produce  a retinal  lesion  In  the  rabbit  fundus  directly 


with  thermocouples.  Their  measurements  of  lesion  temperature  for  a 


10  second  exposure  In  the  rabbit  fundus  Indicated  that  a tem- 
perature Increase  of  15.4°C  above  the  baseline  temperature  of 
37°C  was  required  to  produce  a minimum  visible  lesion. 


They  found  that  the  rise  In  fundus  temperature  required  to 
produce  a visible  lesion  was  approximately  linear  with  pre- 
exposure fundus  temperature  over  a range  from  7°C  below  to 
7.5°C  above  the  normal  37°C  baseline  value. 

The  general  physiological  state  of  the  animal  and 
level  of  anesthesia  can  also  affect  the  retinal  irradiation 
threshold  to  Indeterminate  degrees.  The  condition  of  the  retinal 
vasculature  and  blood  pressure  affect  the  flow  of  blood  in  the 
choriocapi Haris.  Geeraets  and  Ridgeway  (19)  found  a lower 
threshold  value  for  freshly  killed  rabbits  than  for  live 
animals.  Anesthetics  such  as  sodium  pentobarbitol  which  are 
commonly  used  In  threshold  experiments  cause  a partial  decoupling 
of  the  animal's  normal  thermoregulatory  mechanism  and  require 
external  maintenance  of  the  core  temperature  by  heating  pads 
or  heat  lamps. 

Subthreshold  changes  have  been  observed  histologically, 
hlstochemlcally,  and  In  optic  tract  electrical  activity.  McNeer 
et  al.  (20)  found  that  the  ERG  b-wave  was  reduced  In  amplitude 
for  exposures  at  50%  of  threshold  If  sufficient  areas  (approxi- 
mately 44mm2)  of  the  retina  were  Involved.  Other  Investigators 
have  noted  an  Increase  In  the  ERG  for  near  threshold  exposures 
and  a change  In  the  c-wave  at  70%  of  the  threshold  exposure  (21, 
22).  Hempel  (22)  noted  changes  In  the  LGN  and  cortical  responses 
for  near  threshold  values.  Mautner  (23)  observed  a decrease  In 


the  visually  evoked  cortical  potential  at  25%  of  the  ophthal- 
moscoplcally  visible  lesion  threshold.  Farrar  (24)  noted  no 
decrease  In  visual  acuity  for  subthreshold  Injury  to  trained 
monkeys  but  showed  a transient  decrease  In  acuity  for  super- 
threshold lesions. 

In  histological  and  hlstochemlcal  studies,  changes  In 
the  histology,  enzyme  Inactivation  and  protein  denaturatlon 
have  been  observed  for  subthreshold  exposures  to  Intense  light 
sources.  Using  tissue  electrophoresis,  Geeraets  et  al.  (25) 
demonstrated  coagulation  of  retinal  proteins  at  subthreshold 
energies.  They  also  demonstrated  a reduction  in  enzyme  activity 
(DPN  diaphorase)  of  enzymes  which  are  important  in  the  regenera- 
tion of  retinal  photopigments.  DPN  diaphorase,  found  mostly  in 
the  ellipsoids  of  the  visual  receptor  cells,  was  inactivated  as 
evidenced  by  loss  of  staining  ability  for  exposure  energies  10 
to  15  percent  below  threshold. 

Threshold  powers  for  He-Ne  and  A r sources  for  various 
Image  sizes  and  exposure  times  have  been  determined  by  various 
Investigators.  The  data  define  several  Important  variables  which 
Influence  the  response  Including  experimental  animal,  image  size, 
source  wavelength  and  time  of  exposure.  The  monkey  thresholds 
determined  by  Davis  and  Mautner  (14)  are  higher  than  thresholds 
determined  for  the  rabbit.  The  exposure  power  required  to  produce 
a threshold  lesion  Increases  with  decreasing  time  of  exposure. 
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The  total  power  to  produce  a lesion  Increases  with  Increasing 
Image  size  but  not  In  direct  proportion  to  the  area  of  the 
retinal  Image.  As  an  example,  He-Ne  threshold  power  (1)  for  a 
monkey  at  500  milliseconds  exposure  time  Increases  by  a factor 
of  1.67  for  an  Increase  In  the  ratio  of  the  radii  of  1.55. 

Since  the  exposure  area  increase  Is  proportional  to  the  square 
of  the  Increase  in  the  radius,  the  energy  per  unit  area  in  the 
retinal  image  Is  less  for  a threshold  lesion  with  a large  retinal 
image  than  for  a small  Image. 

Temperature  Measurements 

Several  investigators  have  attempted  intraocular  and 
retinal  temperature  measurements  with  thermocouples  and  thermistors 
(26-31).  The  data  from  these  investigations  are  tabulated  in 
Table  I.  Several  experimental  factors  have  caused  the  measurements 
to  be  questioned.  The  measurements  were  made  with  bulk  wire 
thermocouples  of  relatively  large  size  compared  to  the  retinal 
Image  diameter.  Some  of  the  measurements  were  made  during  ruby 
laser  and  other  short  duration  exposure  sources.  The  thermocouples 
In  at  least  one  case  (30)  had  a rise  time  of  40  milliseconds  and 
were  used  to  measure  temperature  rise  for  a 500ysecond  exposure. 

No  analytical  models  were  employed  to  lend  validity  to  the 
accuracy  of  experimental  measurements. 
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Cain  and  Welch  (32)  developed  a technique  for  tempera- 
ture measurement  In  the  rabbit  fundus  which  enabled  accurate 
measurements  of  temperature  transients  in  the  eye.  They  did  not 
measure  the  temperature  associated  with  lesion  thresholds  but 
they  demonstrated  the  validity  of  this  technique.  The  temperature 
probe  consisted  of  a thln-film  copper-nickel  thermocouple  on  the 
tip  of  a quartz  probe-type  substrate.  The  thermal  properties 
of  the  quartz  substrate  approximately  matched  the  properties  of 
the  blcmedium.  The  small  volume  of  the  metal  films  (=  .14%  of 
mass  of  a 20pm  probe  with  .2pm  films  of  copper  and  nickel)  with 
respect  to  the  volume  of  the  quartz  substrate  provided  a sensor 
w 1th  the  approximate  thermal  characteristics  of  quartz . 

The  conductivity  of  fused  quartz  is  .0033  Cal /cm  sec.  °C 
versus  .0015  for  water.  The  densities  of  quartz  and  water  are 
2.29  and  1.0  gm/cm3,  respectively,  and  the  volumetric  specific 
heats  are  0.4  and  1.0  Cal/gmoc.  The  kpc  product  for  quartz  and 
water  are  .00131  Calfycm4  sec.°C2  and  .0015  respectively  (45). 

The  approximate  match  of  the  properties  of  the  sensor  and  the 
biomedium  allows  accurate  transient  temperature  measurements 
for  exposure  durations  longer  than  a few  milliseconds. 


Temperature  Models 


An  evaluation  of  hazards  to  the  eye  presented  by  a 
variety  of  intense  light  sources  requires  an  accurate  thermal  model 
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which  accounts  for  spectral,  temporal  and  spatial  distributions  of 
the  source.  The  development  of  the  model  depends  upon  an  accurate 
knowledge  of  the  spectral  dependence  of  the  transmission  and  ab- 
sorption of  light  in  the  various  layers  of  the  eye  and  the  spatial 
distribution  of  the  Image  in  the  eye.  The  transmission  parameters 
will  be  discussed  in  a later  section. 

Several  models  have  been  developed  which  predict  the 
transient  and  steady  state  temperatures  in  the  eye  (33-41).  The 
models  are  based  on  solutions  of  the  heat  conduction  equation: 


III 

a 3 t 


where  k * thermal  conductivity  (Cal /cm.  sec.°C) 

a * — = thermal  diffusivity  (cm^/sec.) 
pc 

pc  * volumetric  specific  heat  (Cal/cnr.°C) 

T = temperature  rise  (°C) 

A » heat  source  term  (Cal /cm3  sec.) 
p ■ density  (gm/cm3) 
c ■ specific  heat  (Cal/jrPC) 

The  major  differences  among  previously  developed  models 
are  the  assumptions  of  boundary  conditions  and  absorption  phenomena 
In  the  eye.  Vos  (33,  38)  uses  a square,  uniformly  Irradiated 
Image  with  uniform  absorption  In  the  depth  and  width  of  the 
irradiated  volume.  Clarke  et  al.  (35)  have  developed  a steady- 
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0 < e < 2tt 
0 < r < b 
= 0 elsewhere 

where  AQ  is  the  magnitude  of  the  source  term  at  x = 0,  where 

C is  the  depth  of  the  absorbing  layer  and  b is  the  radius 
of  the  disk  Image.  He  assumes  values  of  314  cm-1  and  45  cm-1 

for  the  P.E.  and  choroid  absorption  coefficients  (a)  in  his 
cyllndrlcally  symmetric  two  layer  model.  Calculations  performed 
with  this  model  indicate  a temperature  rise  of  10°C  for  a thres- 
hold  lesion  using  data  of  Ham  et  al.  (5)  and  about  23°C  for  a 
threshold  lesion  from  the  data  of  Allen  et  al.  (42)  for  a white 
light  source. 

Other  steady-state  models  proposed  by  Hansen  et  al. 

(34)  and  Peacock  (43)  are  of  questionable  value  because  they  do 
not  account  for  the  structure  of  the  retina.  They  assume  100% 
absorption  of  the  Input  light  In  a small  thickness  of  the  retina. 

Presently  the  most  comprehensive  model  available  Is  a 
finite  difference  solution  to  the  heat  conduction  equation 
developed  by  Malnster,  White  and  Tips  (39,  40,  41).  The  model 
Includes  the  effects  of  temporally,  spectrally,  and  spatially 
varying  source  terms,  assumes  two  layer  absorption  In  the  P.E. 
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and  choroid,  and  allows  different  conductivities  and  volumetric 
specific  heats  to  be  assigned  to  the  various  media  of  the  eye. 
The  model  calculates  temperature  as  a function  of  position  and 
time  throughout  the  fundus  for  a circularly  symmetric  retinal 


image. 


The  irradiance  (H)  in  the  P.E.  and  choroid  for  a uniform 


image  is  illustrated  in  Figure  2.  The  source  strength  (A)  is  also 
shown  in  the  figure.  Lambert-Beers  (34)  absorption  is  assumed  and 


the  retinal  irradiance  has  the  form, 

H1  a h(r)  HQ  expC-ajZ]  (P.E.) 

H2  s h(r)  HQ  exp(d-)  a2-d]  a^^z)  (Choroid) 

d]  < Z < dj  + d2  where 

h(r)  Is  the  normalized  Irradiance  profile  at  the  retina  and  dj  and 
d2  are  the  depths  of  the  P.E.  and  choroid  respectively. 

The  source  strength  in  each  layer  is  related  to  the 


irradiance  by  the  relationship: 
A. 


1 = 1,2. 


The  various  layers  are  assumed  to  be  homogeneous  and  isotropic  with 
constant  absorption  parameters  and  thermal  parameters  k,  p,  and  c. 
The  retinal  Image  distribution  may  be  entered  directly  into  the 
program  or  calculated  using  the  point  spread  function. 
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Cain  and  Welch  (32)  used  this  model  to  compare  calculated 
and  measured  temperatures.  Measured  retinal  distributions  were 
used  In  the  model  rather  than  calculated  source  characteristics 
and  the  point  spread  function  of  the  eye.  Measurements  Indicated 
that  the  model  predicted  temperatures  were  generally  30%  low 
for  Images  greater  than  10  times  the  size  of  the  thermocouple 
at  exposure  times  greater  than  one  second. 

The  possible  reasons  for  the  differences  between  cal- 
culated and  measured  temperatures  were: 

1)  assumed  values  of  thermal  parameters  k,  p,  and  c. 

2)  assumed  values  for  the  transmittance  of  the  pre- 
retlnal  ocular  media  and  the  absorption  parameters 
for  the  P.E.  and  choroid. 

3)  measurement  error  due  to  direct  absorption  in  the 
thermocouple. 

4)  errors  in  the  measurement  of  the  beam  profile  at  the 
retina. 

5)  errors  In  measurement  of  Incident  corneal  power. 

6)  size  of  the  probe  relative  to  the  size  of  the  image. 

7)  measurement  of  retinal  Image  diameter. 

Of  the  probable  sources  of  error,  Cain  and  Welch  (32)  demonstrated 
that  the  above  points  3 through  6 had  little  effect  upon  the 
measurements.  The  substitution  of  k « .0013  (brain  tissue) 
versus  k ■ .0015  Cal/cm-sec°C  (water)  Increased  the  calculated 
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temperature  for  exposure  times  which  resulted  in  near  steady  state 
temperature  Increases. 

Parameters  in  the  Thermal  Model 

The  model  for  retinal  temperature  rise  which  results 
from  exposure  of  the  ocular  fundus  to  Intense  light  sources 
requires  the  determination  of  thermal  conductivities,  volumetric 
specific  heats  and  light  absorption  in  the  various  layers  of  the 
eye.  Generally,  the  values  for  conductivity  (k)  and  volumetric 
specific  heat  (pc)  are  taken  to  be  the  same  as  for  water  in  each 
layer  in  the  model.  Only  the  absorption  of  light  in  the  different 
layers  of  the  eye  is  variable.  However,  a major  source  of  error 
in  the  calculation  of  temperature  for  long  exposures  is  the  value 
assumed  for  conductivity  and  the  volume  of  tissue  heated  which  is 
related  to  the  absorption  coefficients  of  the  layers  within  the 
retina.  For  short  duration  exposures  and  for  tissue  rise  time 
calculations,  the  volumetric  specific  heat,  the  absorption  coef- 
ficients and  the  measurement  of  retinal  irradlance  affect  the 
model  solutions. 

The  conductivity  value  for  water  is  1.5  x 10"  3 cal  cm"1 
sec”1  °C_1  for  human  blood  1.21  x 10”3  cal  cm'1 sec”1  °C-1  and  for 
rat  blood  1.26  x 10”3  cal  cm”1  sec”1  °C-1  (44).  From  the  conduction 
equation  at  steady  state: 
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which  indicates  that  the  steady  state  temperature  is  inversely 
proportional  to  the  conductivity  and  directly  proportional  to  the 
volume  of  tissue  heated  by  the  source  term  A.  From  the  values  of 
conductivity  given  above,  the  values  for  blood  would  imply  a 
higher  steady  state  temperature  than  the  value  for  water. 

The  source  term  (A)  includes  an  absorption  coefficient 
a and  is  given  by: 

ApE  = h(r)  H(^1  exp  [-a-j z]  0 <z  < d1 

in  the  P.E.  where  d^  equals  the  depth  of  the  P.E.  and  ^ equals 
the  absorption  coefficient  of  the  P.E. 

For  the  choroid,  the  source  term  is  given  by: 

Ach  = h(r)  Hoa2exptdla2  " dlal  -a2zl 
d-j  < z < d2  + d-j 

where  a j and  d£  are  the  absorption  coefficient  and  depth  of  the 
choroid  respectively.  The  volume  of  tissue  heated  enters  Into 
the  effective  rise  time  of  the  temperature-time  history  at  a 
point  In  the  tissue  and  Into  the  final  temperature  at  that 
point. 

For  very  short  times  for  which  no  significant  conduction 
has  taken  place,  the  conduction  equation  is  given  by: 
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and  the  temperature  is  approximately: 


T * t, 
pc 


at  z = 0 and  r = 0. 


So,  at  short  times,  the  slope  of  the  temperature  versus  time  is 
proportional  to  aHQ  and  inversely  proportional  to  pc.  These 
simplifying  assumptions  show  the  importance  of  the  absorption 
coefficients  in  the  calculation  of  the  temperature  versus  time 
profiles  in  the  tissue. 


Measurement  and  Estimation  of  Absorption  Parameters 


The  absorption  and  transmission  of  the  various  media 
within  the  eye  have  been  measured  by  several  different  investiga- 
tors. The  measurements  of  Geeraets  et  al.  (7,  8)  and  Boettner 
and  Wolter  (9)  are  the  most  frequently  referenced  values  in  the 
literature.  The  measurements  of  absorption  and  transmission 
have  been  accomplished  In  vitro  on  excised  human,  primate  and 
rabbit  eyes.  Calculations  of  absorption  parameters  for  the 
assumption  of  Lambert- Beers  absorption  have  been  made  by  several 
different  Investigators.  A summary  of  past  calculations  is  pre- 
sented In  Table  II.  Since  few  separate  measurements  for  the 
pigment  eplthelulm  and  choroid  exist  In  the  literature,  the 
general  assumption  has  been  made  that  50%  of  the  light  reaching 


the  P.E.  Is  absorbed  In  that  layer.  Generally,  the  absorption 
coefficients  of  Table  II  do  not  follow  this  assumption  as 
illustrated  in  Figure  3. 

In  Figure  3 is  Illustrated  a plot  of  the  percent  ab- 
sorption versus  az  for  the  assumption  of  Lambert-Beers  absorption 
in  the  tissue.  For  a 50%  absorption  of  light  in  the  10pm  thick- 
ness of  the  P.E.,  az  Is  .69  this  implies  that  a,  the  absorption 
coefficient  for  the  P.E., is  690  cm~^  or  less.  In  Table  II,  three 
of  the  four  tabulated  values  exceed  this  value.  In  fact,  for  the 
value  of  2080  cm-^  for  a in  the  P.E.,  88%  of  the  light  incident 
on  this  layer  is  absorbed. 


TABLE  II:  Values  of  Absorption  Coefficients 
Calculated  by  Past  Investigators. 


Reference 

Animal 

TOM 

APE 

ACH 

Sclera 

White  (45) 

Monkey 

.79 

2080 

208 

.044 

MCV  (by  White) 
(45) 

Rabbit 

.82 

832 

83.2 

.770 

Boettner  (Brooks 
AFB)  (44) 

Monkey 

.51 

1622 

189 

Tom  White  (45) 

Rabbit 

.636 

637.6 

76.28 

TOM  = The  transmission  of  the  preretlnal 

ocular 

media. 

APE  » The  absorption  coefficient  of  the  P.E. 

ACH  = The  absorption  coefficient  of  the  choroid. 


FIGURE  3:  Percent  absorption  In  an  absorbing  layer 
for  Lambert-Beers  model . The  percent  of 
light  reaching  the  front  of  the  layer  which 
Is  absorbed  In  that  layer  is  plotted  versus 
the  absorption  coefficient  (a)  times  the 
tissue  depth  (z 
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It  Is  also  apparent  that  large  discrepancies  in  the 
transmission  of  the  preretlnal  ocular  media  exist  In  the  litera- 
ture from  the  values  of  TOM  given  In  the  table.  The  transmission 
of  the  preretlnal  ocular  media  Includes  absorption  in  the  tissues 
and  fluids  anterior  to  the  P.E.  and  reflection  from  boundaries 
of  different  refractive  Index.  Major  differences  in  the  values 
determined  for  TOM  exist  between  the  measurements  of  Geeraets  and 
Berry  (7,  8)  and  Boettner  (9).  Geeraets  and  Berry  measured  both 
directly  transmitted  light  and  scattered  light  while  Boettner 
attempted  to  separate  the  direct  and  scattered  components  of 
light.  Boettner' s assumption  was  that  light  within  a 1°  cone 
at  the  exit  from  the  tissue  being  measured  was  directly  transmitted 
and  the  light  outside  this  cone  was  scattered  light.  The  difference 
in  the  two  sets  of  measurements  is  approximately  30%. 

As  light  impinges  upon  the  eye,  it  encounters  first  an 
air-water  interface.  A portion  of  the  light  is  reflected  due  to 
the  different  indices  of  refraction  according  to  the  relationship: 
(n]  - n2 

R = o 

(nj  + n2r 

where  n^  and  n^  are  the  refractive  indices  of  the  two  media  and  R 
is  the  fraction  of  incident  power  reflected  (47).  This  relation 
holds  only  for  multiple  surfaces  between  media  which  are  separated 
by  distances  greater  than  a wavelength.  For  very  thin  films  such 
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as  water  on  the  cornea,  the  actual  reflection  may  differ  markedly 
from  this  value  due  to  Interference  effects.  Also,  the  use  of  a 
contact  lens  presents  a higher  index  than  water  to  the  Incoming 
light.  Contact  lenses  are  frequently  used  in  experiments  to 
prevent  corneal  drying.  From  the  formula  above,  the  reflection 
of  light  from  the  various  surfaces  of  the  eye  may  be  calculated. 

A summary  of  these  effects  is  presented  in  Table  III. 


TABLE  III:  Power  Reflection  at  Various  Ocular  Interfaces 

for  X = 500oX 

Material  or  Index  of  Refraction  Index  Anterior  % Reflection 
Ocular  Medium  of  Medium  (48)  to  Medium 

Glass  (contact 


lens) 

1.5 

1.00  (air) 

4.0 

Cornea  (contact 
lens) 

1.37 

1.5  (contact) 

.2 

Cornea  (air) 

1.37 

1.00  (air) 

2.44 

Aqueous  humor 

1.336 

1.37 

.02 

Crystal ine  lens 

1.4 

1.336 

.06 

Vitreous  humor 

1.336 

1.4 

.06 

Neural  retina 

— 

1.336 

3.0 

TOTAL  REFLECTIONS 
(with  contact) 

7.34 

(no  contact  lens) 

5.58 

- i.  ~ TIm. 


_ T . ut 
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The  reflection  of  the  neural  retina  is  from  Boettner  (9)  at  500o)l 
wavelength.  For  a contact  lens  on  the  cornea,  approximately 
7.3*  of  the  incident  light  is  reflected  at  the  various  surfaces 
of  the  eye.  With  no  contact  lens,  approximately  5.6*  is  reflected. 
Assuming  an  absorption  coefficient  of  an  equal  length  of  water  for 
the  preretinal  ocular  media  ( a»  .05  cnr^  (35)  ),  the  fraction  of 
corneal  light  reaching  the  pigment  epithelium  is  given  by. 


TOM  = 1 - * Reflected  -e 


where  z is  approximately  1.8  cm  for  the  monkey  eye.  Therefore, 

TOM  is  approximately  83-85*. 

As  shown  in  Table  II,  all  measured  values  are  lower  than 
the  calculated  values.  One  reason  is  that  the  eye  is  not  a homo- 
geneous medium  and  cannot  be  approximated  accurately  by  an 
equivalent  length  of  water.  Some  light  is  scattered  out  of  the 
image  due  to  Inhomogeneities  In  the  ocular  media  and  nonlinear 
optical  distortions  in  the  optical  system. 

A comparison  of  the  results  presented  In  Table  II, 
however,  show  striking  differences  In  transmission  values  for  the 
monkey  determined  by  two  groups  of  Investigators.  The  value 
obtained  by  Boettner  for  direct  transmission  is  much 
lower  than  the  value  from  Geeraets  and  Berry.  Boettner 
also  obtained  a total  transmission  figure  of  65-70*  for 
the  monkey  (9).  From  the  above  comparisons,  it  can  be  seen  that 


_ - _|  ■-  — - 


■up  ,v  ■ ■ m 


m 


28 


even  one  of  the  seemingly  simpler  measurements  yields  large  dif- 
ferences between  Investigators.  The  measurement  of  TOM  is  a 
critical  parameter  in  the  model  for  retinal  temperature  rise  since 
the  magnitude  of  the  source  term  vaires  directly  with  changes  in 
TOM.  Since  the  model  is  a linear  equation,  the  temperature  varies 
linearly  with  changes  in  the  value  for  preretinal  ocular  media 
transmission. 

The  absorption  parameters  in  Table  II  also  appear  to  be 
essentially  unknown.  There  are  few  measurements  of  the  P.E.  and 
choroid  separately,  but  those  measurements  indicate  a nearly  50-50 
split  for  rabbits  and  humans  in  the  absorption  of  light  in  the 
P.E.’  and  choroid  (7,  8).  From  Figure  3,  it  is  apparent  that  t(^ 
absorption  coefficients  listed  in  the  table  are  too  high  for  the 
P.E.  since  a 10pm  thickness  even  for  a = 1622  cm'1  results  in  80% 
of  the  light  being  absorbed  in  the  P.E.  From  Geeraets  and  Berry's 
(7,  8)  measurements,  approximately  11%  of  the  light  at  488oJi  is 
transmitted  through  the  P.E.  and  choroid.  With  a 50%  absorption 
of  this  light  in  each  layer,  44.5%  of  the  light  reaching  the  P.E. 
is  absorbed  in  the  P.E.  and  80.2%  of  the  light  reaching  the  choroid 
is  absorbed  in  the  choroid.  From  Figure  3,  az  is  approximately  .60 
for  the  P.E.  and  which  results  in  absorption  coefficients  of  600 
cm-1  and  162  cm"1  for  the  P.E.  and  choroid,  respectively.  The  50% 
split  of  absorbed  radiation  is  somewhat  arbitrary  and  the  choice  of 
TOM  is  also  open  to  question.  The  total  transmission  values  may  be 


.... 


■BHBBnHSaBBBHBOl 


■ a 


29 


too  high  and  the  direct  transmission  values  are  almost  certainly 
too  low.  The  choice  of  TOM  Is  therefore  In  the  range  from  51  to 
80  percent  for  the  monkey  eye.  A choice  of  65  percent  transmission 
should  therefore  be  within  ± 23*  of  the  correct  answer  for  the 
transmission  of  the  preretlnal  ocular  media  and  Is  In  close 
agreement  with  the  total  transmission  value  of  Boettner  (9).  The 
change  In  TOM  directly  affects  the  calculations  of  retinal 
temperature;  a 20*  change  In  TOM  causes  a 20*  change  In  the 
calculated  temperature  Increase. 

Probe  Measurement  of  Retinal  Temperature 

The  Invasive  measurement  of  tissue  temperature  with  a 
probe- type  thermocouple  requires  a reasonable  model  for  the  pertur- 
bation of  the  measurement  by  the  sensor.  The  probe  presents  several 
problems  In  temperature  measurement  some  of  which  may  be  accounted 

9 

for  In  analysis,  others  which  can  only  be  estimated.  In  order  to 
make  accurate  temperature  versus  time  measurements,  the  rise 
time  of  the  probe-tissue  system  must  be  tissue  limited.  This  Im- 
plies that  the  probe  step  response  In  the  medium  must  be  faster 
than  the  rise  time  of  the  medium  for  the  measurement  conditions. 

A common  method  used  to  characterize  the  ability  of  a thermocouple 
sensor  to  measure  dynamic  temperature  changes  Is  the  Idealized 
step  response.  This  response  Is  not  only  a probe  characteristic 
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but  it  also  depends  upon  the  medium  in  which  the  probe  is  imbedded 
as  will  be  demonstrated  later. 

The  thermal  resistance  of  the  probe  to  axial  conduction 
down  its  leads  and  the  thermal  resistance  from  the  probe  to  the 
tissue  are  also  very  important  properties  of  the  probe-tissue 
system.  If  the  axial  conduction  resistance  is  lower  than  or  low 
enough  with  respect  to  the  probe  to  tissue  Interface  resistance, 
the  probe  may  measure  temperatures  which  are  significantly  lower 
than  the  temperature  of  the  medium.  For  small  cross-sectional 
lead  areas  this  conduction  resistance  can  be  maximized  in  a real 
system. 

In  the  case  of  the  measurements  for  this  research,  the 
direct  absorption  by  the  probe  of  Incident  radiation  must  be  small 
compared  to  the  absorption  In  tissue.  The  differential  absorption 
of  radiant  energy  by  the  probe  and  tissue  gives  rise  to  a direct 
absorption  component  of  temperature  rise  which  tends  to  make  the 
probe  indicate  a temperature  which  Is  higher  than  the  temperature 
of  the  surrounding  tissue.  Since  the  rise  time  of  the  direct 
absorption  component  is  much  faster  than  convection  from  tissue, 
the  measured  rise  time  may  be  separated  Into  two  components. 

The  rise  time  or  time  to  nearly  steady-state  temperature 
depends  upon  the  volume  of  tissue  heated.  A smaller  volume  will 
reach  stea<ty-state  faster  than  a larger  volume.  For  the  reasons 
stated  above,  the  volume  of  tissue  heated  should  be  much  larger 
than  the  volume  occupied  by  the  probe.  Since  the  tissue  in  this 
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research  Is  heated  by  a radiant  light  source,  the  shadowing  by 
the  probe  will  remove  less  than  one  percent  of  the  source  term  If 
the  radius  of  the  Image  Is  10  times  that  radius  of  the  probe. 

Axial  conduction  down  the  probe  will  also  be  minimized  for  larger 
Images. 

The  characteristics  of  the  probe  to  a step  Increase  In 
temperature  have  been  discussed  by  Reed  (49)  and  Cain  and  Welch  (32). 
Reed  calculated  the  rise  time  of  a cylinder  In  an  Infinite  medium  for 
various  simplifying  assumptions.  First,  the  probe  was  entirely 


embedded  In  the  medium  with  no  axial  flow.  Second,  the  heat  trans- 
fer between  the  probe  and  medium  was  by  conduction.  The  problem 
which  Reed  solved  was  for  a probe  Initially  at  unit  temperature 
and  at  time  t 3 0 It  was  placed  In  an  Infinite  medium  at  zero 
temperature.  This  solution  can  also  be  found  In  Carslaw  and  Jaeger 
(50).  The  calculated  fall  time  (cooling  time)  Is  equivalent  to 
heating  a probe  In  a medium  at  elevated  temperature. 

Reed  numerically  evaluated  the  solution  from  Carslaw  and 
Jaeger  (50)  and  obtained  plots  of  non-dimensional  rise  time  t versus 
relative  volumetric  specific  heats.  The  results  of  his  analysis 
are  tabulated  for  a few  materials  of  Interest  to  this  research  In 
Table  IV.  The  dimensionless  rise  time  Is  given  by: 


where; 

a = k/pc  * thermal  dlffuslvlty, 

A = radius  of  the  thermocouple  probe, 
c 

and,  t = real  time. 

For  the  case  of  a probe  In  an  Infinite  medium,  the  dlffuslvlty  Is 
the  dlffuslvlty  of  the  medium.  The  relative  volumetric  specific 
heat  used  In  the  determination  of  rise  time  is: 

* = pc  medium 
pc  probe 

where  pc  of  the  probe  Is  assumed  to  be  the  value  for  quartz. 

TABLE  IV:  Relative  Volumetric  Specific  Heats,  Rise 
Times  and  Diffusivities  for  Several  Media, 
from  Reed  (49) 


Medium 

v T1 0-90% 

o 

a (cm  /sec) 

Water 

2.43  .47 

1.44  x 10-3 

Air 

8.8  x 10'4  6 x 102 

.187 

Biomedium 

-A 

(skin) (45) 

2.22  .64 

7.4  x 10  * 

In  Figure  4 Is  Illustrated  the  rise  time  (10-90*)  response  of  quartz 
cylinders  embedded  In  the  media  listed  in  Table  IV.  Experimental 
points  from  Cain  and  Welch  (32)  are  also  illustrated  In  the  Figure. 
The  measurements  were  made  by  Cain  by  driving  thermocouples  Into  a 


heated  water  bath.  The  measured  values  Indicate  a slightly 
longer  rise  time  than  the  calculated  values  for  all  probe  sizes. 

This  Is  to  be  expected  since  the  model  assumes  Intimate  thermal 
contact  between  the  probe  and  the  medium.  The  model  Indicates  a 
rise  time  (10-90%)  for  the  10pm  radius  probe  in  water  of 
3.26  x 10"*  seconds  while  the  experimental  value  is  4.8  x 10"* 
seconds.  For  a 30ym  radius  probe  in  water  the  model  and  experiment- 
al values  are  2.93  x 10'^  and  3.6  x 10"^  seconds,  respectively. 

The  model  predicts  a rise  time  in  air  which  is  approximately 
100  times  the  rise  time  In  water  and  37  times  the  rise  time  in 
the  biomedium. 

Probe  to  tissue  thermal  resistance  was  assumed  to  be 
zero  in  the  above  analysis.  Reed  (49)  also  performed  the  analysis 
for  cases  for  which  thermal  resistance  between  the  probe  and 
tissue  was  present.  The  surface  resistance  to  conduction  y was 
used  in  this  analysis.  This  resistance  is  defined  by: 
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y = kR/£c  (cm"2) 
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0* 


where  k Is  the  thermal  conductivity,  R Is  the  surface  resistance 
per  unit  of  surface  area  of  the  sensor  and  lc  is  a characteristic 
dimension  of  the  sensor  (diameter  or  length).  For  a "typical" 
value  of  a less  than  ideal  thermal  Interface  (y  * 4),  the  response 
time  increases  by  a factor  of  4.4  for  quartz  in  a biomedium  and 
1.8  for  quartz  in  air.  Thus,  for  a realistic  situation,  the 
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response  time  of  the  model  for  quartz  in  a biomedium  for  a 10pm 
radius  probe  is  3.8  milliseconds.  One  would  expect  from  Reed's 
analysis  and  Cain's  measurements  that  for  times  greater  than 
10  milliseconds,  the  probe  response  is  essentially  complete 
and  the  perturbation  of  temperature  measurement  due  to  probe 
transient  response  is  negligible. 

Rate  Process  Models 

The  biological  damage  which  results  from  elevated 
temperature  in  tissue  has  been  modeled  by  several  Investigators. 
Henrlques  (51)  was  the  first  to  apply  a rate  process  model  to 
thermal  damage  to  a biological  tissue,  porcine  epidermis. 

Fugltt  (52)  and  Stoll  (53,  54,  55)  modified  the  model  with 
empirically  derived  rate  constants.  Hu  and  Barnes  (56),  Kach  (57), 


A Is  a pre-exponential  or  frequency  factor,  EQ  Is  the  experimental- 
ly derived  activation  energy,  T Is  the  absolute  temperature  (°K) 
and  R Is  the  gas  constant.  The  reaction  velocity  constant  may  be 
written  as  a function  of  a change  In  free  energy,  AF  where: 


AF  * AH  - T A S 


In  which  AH  Is  the  barrier  heat  energy  and  AS  Is  the  barrier 
entropy.  In  the  equation  for  AF,  the  temperature  change  AT  is 
assumed  to  be  zero.  The  reaction  rate  constant  thus  becomes: 

k,  _ kT  -AF/RT  kT  AS/R  -AH/RT 
K ~ - — e = - — e » 


where  k = Boltzmann's  constant,  and  h * Planck's  constant. 
The  derivative  of  In  k'  with  respect  to  T is: 


d (In  k')  = AH  + RT 

ciT  7 

RTZ 

which  is  Identical  to  the  Arrhenius  equation  with  EQ  * AH  + RT. 

The  reaction  kinetics  for  a first-order  reaction  give  for 
the  change  in  concentration  with  respect  to  time: 


dc  . , 
3t  “ k c 


where  c Is  the  concentration.  The  fraction  of  remaining  molecules 
at  time  t.  Is  given  by. 
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7 ~ * e_/o  k*  dt 
co 


The  surviving  fraction  Is  therefore  dependent  upon  the 
barrier  heat  energy  of  the  reaction,  the  temperature  and  the 
length  of  time  at  an  elevated  temperature. 

The  term: 

kT  A +AS/R 
TT  e 

Is  assumed  to  be  constant  with  temperature,  even  though  it  does 
vary  slightly.  This  assumption  is  not  too  bad  since  AS  is 
relatively  Insensitive  to  temperature  and  for  typical  proteins 
AS/R  is  * 30  to  250  (55). 

If  k*  Is  relatively  independent  of  temperature,  a plot 
of  In  S versus  time  yields  a line  with  slope  -k',  since 


In  (c/c  ) = -k'  t 
o 

Is  of  the  form  y * mx  + b.  Plotting  AF  versus  T allows  determina- 
tion of  AH  and  AS  since: 


1F  ■ -RT  *"  <ir7h> 


AF  « AH  - T AS 


s > act  »*.  • rnm  -i 
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if  the  variation  of  k'  with  temperature  is  known.  Thus,  the 
various  constants  can  be  determined  from  a plot  of  experimentally 
determined  data. 

The  applications  of  a rate  process  models  to  biological 
tissues  have  been  based  on  rather  arbitrary  criteria.  No  quanti- 
tative assays  of  protein  or  enzyme  deactivation  associated  with 
threshold  burns  have  been  made.  Generally,  two  criteria  are 
assumed,  the  first  is  complete  cellular  necrosis  associated  with 
a "damage"  integral  (51,  52): 

ft  - /*  A.  "E/RT  dt 
o v® 

where  ft  * 1 at  threshold. 

The  second  commonly  used  criterion  is  for  S,  the  surviv- 
ing fraction  equal  to  e“l  which  corresponds  to: 

c ■ c0/e  = 368  V 

Typical  values  of  the  rate  constants  temperatures  and 
activation  energies  and  barrier  entropies  for  denaturation  of 
proteins  and  enzyme  Inactivation  have  been  documented  by  Wood  (59). 
The  value  for  AH  Is  between  35,600  and  198,000  for  the  enzymes  and 
proteins.  The  temperatures  of  denaturation  and  Inactivation  vary 
from  25°C  to  120°C. 


Significance 


Many  questions  must  be  answered  if  a comprehensive 
theory  of  retinal  Injury  Is  to  be  developed.  Ideally,  a single 
model  could  be  used  to  predict  Injury  to  ocular  tissues  for  a 
variety  of  exposure  conditions  including  wavelength,  pulse 
width,  pulse  repetition  frequency  and  image  geometry.  This  may 
not  be  a practical  goal  with  respect  to  a simple  thermal  model, 
but  may  need  to  Include  the  short  and  long-term  effects  mentioned 
earlier.  The  thermal  model  may  be  of  some  practical  use,  however, 
when  coupled  with  a rate  process  model  over  a limited  range  of 
exposure  conditions  (say  100pm  seconds  to  10  seconds).  It  may  be 
of  use  in  clinical  applications  for  the  prediction  of  the  extent  of 
damage  In  retinal  diathermy.  Also,  the  development  of  safety 
standards  could  be  aided  by  the  thermal  model  coupled  with  actual 
temperature  measurements  in  the  retina  or  other  ocular  tissues. 


CHAPTER  III 
EXPERIMENTAL  PROCEDURE 


Surgical  Procedure 

A surgical  procedure  has  been  developed  which  exposes 
the  backside  of  the  eye  of  anesthesized  monkeys  (60).  Nembutal 
is  injected  through  a catheter  inserted  in  a large  saphenous  vein 
to  keep  the  animal  sedated  throughout  the  surgery  and  the  ex- 
periment. Earbars  mounted  on  a platform  secure  the  animal's  head. 

A midline  incision  in  the  scalp  is  made  from  the  supra- 
orbital margin  of  the  frontal  bone  extending  to  approximately 
3 cm  beyond  the  junction  of  the  frontal  and  parietal  bones. 

Two  incisions  in  the  skin  are  made  approximately  perpendicular 
to  the  midline  incision.  One  along  the  supra-orbital  margin 
from  the  midline  to  the  zygomatico-temporal  foramen.  The 
other  incision  extends  from  the  midline  incision  to  the  back  of 
the  ear.  The  fascia  of  the  temporal  muscle  is  cut  along  the 
superior  temporal  line  and  the  temporal  muscle  is  elevated  and 
reflected  away  from  the  skull  to  expose  the  posterior  portion  of 
the  bony  orbit.  A dental  drill  Is  used  to  remove  bone  to  expose 
the  back  of  the  eye  In  an  area  between  the  anterior  and  posterior 
ethmoidal  foramena.  The  perl -orbital  fat  and  a portion  of  the 
lacrimal  gland  are  then  removed  and  the  sclera  Is  cleared  of  all 
fascia  for  probe  Insertion.  The  insertion  area  for  the  probe  is 
generally  within  1-2  millimeters  of  the  macula  or  in  the  macula. 
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Three  sutures  are  placed  In  the  conjuctlva  to  bind  the 
eye  to  an  eyeholder  (a  ring  connected  to  an  adjustable  arm 
mounted  on  the  animal  platform)  which  prevents  eye  movement  during 
temperature  measurements.  The  eyeholder  Is  adjusted  to  the  correct 
height  and  position,  and  a contact  lens  Is  mounted  on  the  cornea 
to  prevent  drying  of  the  cornea. 

The  animal  platform  Is  positioned  with  a five  degree  of 
freedom  animal  holder  to  provide  a clear  view  of  the  retina  through 
a fundus  camera.  A point  as  near  as  possible  to  the  macula  Is 
chosen  as  the  site  of  probe  Insertion.  This  point  is  positioned 
In  the  crosshairs  of  the  fundus  camera,  and  a low  intensity  Argon 
C.W.  laser  beam  (4880?)  Is  directed  to  It.  The  size  and  intensity 
of  the  laser  Image  at  the  retina  Is  controlled  with  lenses  and 
neutral  density  filters.  The  resulting  laser  image  can  be  observed 
through  the  sclera  and  serves  as  a marker  for  insertion  of  the 
probe.  Low  laser  energies  are  used  during  probe  insertion  to 
prevent  thermal  damage  to  the  fundus. 

Thermocouple  Probe  Fabrication 

The  thermocouple  probes  used  In  this  study  were  developed 
In  the  Engineering  Mechanics  Laboratory  at  the  University  of  Texas 
by  Professor  E.A.  Rlpperger  and  Mr.  W.D.  Martins.  The  probes  were 
manufactured  for  this  study  by  W.D.  Martins.  The  probe  fabrication 
procedure  and  characteristics  of  the  probe  have  been  described 
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elsewhere  (61,  62).  Only  a brief  description  of  the  probe  geo- 
metry and  characteristics  will  be  presented  In  this  section. 

The  thermocouple  probe  is  a copper  nickel  junction 
vacuum  deposited  on  a quartz  substrate.  The  quartz  substrate  is 
drawn  to  a tip  diameter  of  10-20ym  from  a millimeter  quartz  rod. 
After  the  substrate  is  cleaned,  the  nickel  lead  is  attached  to 
the  deposited  film.  The  tip  of  the  probe  is  masked  and  an 
Insulating  material  Is  deposited  on  the  probe  and  lead  with 
the  exception  of  the  probe  tip.  A copper  film  is  then  deposited 
over  the  probe  and  tip,  forming  a junction  only  at  the  tip  of  the 
probe.  A copper  lead  Is  attached  to  the  film  and  a final  coat 
of  Insulator  Is  deposited  over  the  entire  probe,  including  the 
tip,  to  prevent  Interaction  of  the  metal  junction  with  body  fluids. 

The  copper  nickel  thermocouple  junction  provides  an  EMF 
(electromotive  force)  of  18-21yV  per  degree  Celsius.  The  probes 
are  linear  in  response  to  water  bath  temperatures  over  the  range 
from  0°C  to  100°C.  Typical  rise  times  (10-90X)  of  a lOym  probe 
In  a water  bath  are  less  than  1 millisecond.  The  probe  on  the 
quartz  substrate  Is  mechanically  strong  enough  to  allow  insertion 
Into  soft  tissue  media  without  breaking. 

Instrumentation 

A block  diagram  of  the  experimental  equipment  Is 
Illustrated  In  Figure  5.  The  thermocouple  probe  Is  attached  to  a 


FIGURE  5:  Block  diagram  of  recording  system 


holder  which  Is  held  In  position  by  a micromanipulator  which 
allows  three  axis  movement.  Axial  motion  of  the  probe  is  accom- 
plished with  a hydraulic  drive  which  allows  2pm  or  larger 
movement.  The  probe  Is  attached  to  an  Instrumentation  amplifier 
and  a reference  junction  with  copper  and  nickel  leads.  The 
amplifier  and  the  reference  junction  are  maintained  at 
37°C  1 .05°C  in  a proportionally  controlled  component  oven. 

The  thermocouple  preamplifier  is  nulled  to  37°C  by  placing  the 
thermocouple  in  a 37°C  bath.  Small  temperature  changes  from  the 
reference  (body)  temperature  may  thus  be  measured.  The  gain  of 
the  preamplifier  is  1000  and  its  bandwith  is  DC  to  65KHz.  The 
amplifier  Input  noise  level  for  a typical  probe  is  4yV  RMS. 

The  output  of  the  thermocouple  preamplifier  is  fed  to  a 
Brush  Universal  Amplifier  with  DC  to  lOKHz  bandwidth.  This  ampli- 
fier provides  variable  gain  for  the  signal  with  ±1%  linearity. 

The  output  of  the  Brush  amplifier  is  recorded  on  analog  tape  for 
digitizing  and  on  a Clevlte  Brush  Strip  Chart  recorder  to  provide 
real  time  Information.  The  analog  tape  recorder  has  a bandwith 
of  DC  to  14KHz  at  the  recording  speed  of  30  Inches  per  second. 

Experimental  time  control  and  laser  shutter  operation  is 
provided  by  a digital  timer.  The  timer  provides  digitizing  pulses 
and  timing  control  to  vary  the  exposure  duration  by  opening  and 
closing  a mechanical  shutter.  The  shutter  provides  a variable 
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duration  exposure  with  approximately  .2  millisecond  opening  time. 
Digitizing  pulses  and  shutter  pulses  are  recorded  on  the  analog 
tape  with  the  temperature  Information  for  later  data  reduction. 

An  EG&G  Radiometer  Is  used  to  monitor  the  laser  power 
for  each  exposure.  The  radiometer  Is  calibrated  by  an  NBS  trace- 
able thermopile  with  an  accuracy  of  ±5%.  A light  pulse  signal 
from  the  radiometer  Is  recorded  for  timing  reference  on  the 
analog  tape. 

A parallel  real  time  monitor  of  the  temperature  rise 
Is  provided  with  a 1 MHz  bandwlth  oscilloscope.  Since  the  chart 
recorder  bandwidth  Is  limited  to  90Hz,  accurate  real  time  deter- 
mination of  the  direct  absorption  of  the  thermocouple  may  be 
accomplished  by  the  scope  monitor  and  polarold  photographs.  The 
scope  provides  a real  time  feedback  for  positioning  the  probe  in 
the  peak  of  the  temperature  profile. 

Visual  observation  of  lesion  formation  is  accomplished 
with  a fundus  camera  mounted  on  the  experimental  table.  A reticle 
In  the  fundus  camera  eyepiece  Is  calibrated  to  allow  visual  measure 
ment  of  the  lesion  size.  The  appearance  of  a lesion  five  minutes 
post  exposure  Is  the  criterion  for  threshold.  Fundus  photographs 
of  lesions  and  probe  Insertion  sites  may  also  be  taken  with  the 
aid  of  the  fundus  camera. 


W.'ilf  - 


•*.  vmm  i 


46 


1 1 1 


Temperature  Measurements 

A micro- thermocouple  probe  (32)  Is  Inserted  Into  the 
sclera  by  the  use  of  an  ultrasonic  dentral  drill.  A hole  Is 
first  drilled  in  the  sclera  with  a sharply  pointed  tungsten 
probe  (l-3ym  tip  diameter).  Then  a long  focal  length  microscope 
mounted  above  the  animal  Is  used  to  locate  the  hole  for  the  In- 
sertion of  the  10  to  20um  temperature  probes.  After  Insertion  of 
the  probe.  Its  depth  is  adjusted  by  observing  the  position  of  the 
probe  tip  In  the  vitreous  humor  with  the  fundus  camera.  When 
the  probe  is  just  visible.  It  Is  approximately  10-lOOym  in  front 
of  the  plgnent  epithelium.  With  the  probe  visible,  the  animal 
holder  Is  adjusted  until  the  probe  tip  Is  In  the  approximate 
center  of  the  retinal  Image.  The  direct  absorption  characteristic 
of  the  thermocouple  during  1.6  millisecond  shuttered  laser  pulses 
Is  used  to  position  the  probe  In  the  center  of  the  beam  (point  of 
maximum  temperature  Increase  due  to  direct  absorption)  (60).  By 

. 

rotating  the  animal  about  the  vertical  axis  of  the  eye  (a  hori- 
zontal scan),  the  relative  retinal  Image  Irradlance  profile  and 
image  radius  are  measured  by  the  direct  absorption  of  laser  energy 
by  the  thermocouple.  This  measurement  Is  readily  quantified 
because  the  temperature  response  time  due  to  direct  absorption  Is 
much  faster  than  that  due  to  conducted  heat  (32).  Temperature 
changes  In  the  probe  resulting  from  direct  absorption  of  Incident 
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radiation  appear  as  a step  Increase  with  the  onset  of  the  light 
pulse  and  a step  decrease  at  the  end  of  the  light  pulse. 

The  probe  Is  withdrawn  to  the  center  of  the  pigment 
epithelium  from  the  vitreous  In  2-1 Opm  steps.  The  position  of 
maximum  temperature  rise  due  to  conduction  for  Irradiations  of 
ten  milliseconds  Is  theoretically  In  or  near  the  pigment  epithelium. 
From  this  reference  point  In  the  center  of  the  laser  beam  and  the 
axial  position  of  highest  temperature  rise  measurements  can  be 
made  radially  and  axially  In  the  fundus.  The  position  of  the 
temperature  probe  with  respect  to  the  retinal  Image  may  be  varied 
In  the  horizontal  plane  of  the  eye  In  4pm  to  100pm  steps  to  provide 
radial  measurements  (horizontal  scan).  Measurement  may  be  made 
along  the  path  of  the  probe  of  Insertion  In  2pm  or  larger  steps 
to  provide  temperature  measurement  In  the  axial  direction. 

Generally*  the  probe  track  is  at  an  angle  of  10°  to  30°  for  macular 
Insertions.  Fundus  temperatures  are  measured  for  variations  In  light 
duration  from  .0016  second  to  10  seconds.  All  measurements  repre- 
sent temperature  variation  from  a reference  value  of  37°C. 

Retinal  temperature  Is  monitored  and  maintained  at  37°C 
± 1°C  for  all  threshold  measurements.  The  37°C  baseline  Is  set 
prior  to  probe  Insertion  by  Inserting  the  probe  In  a water  bath 
maintained  at  37°C.  The  baseline  Is  set  on  a Brush  strlpchart 
recorder  and  the  null  adjustment  on  the  thermocouple  amplifier  Is 
set  to  give  0 volts  output  for  the  probe  In  the  37°C  bath. 
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An  Infrared  lamp  is  used  to  maintain  the  fundus  at  37°C  throughout 
the  experiment.  The  sclera  Is  kept  moist  with  cotton  soaked  In 
saline  which  Is  positioned  around  the  thermocouple  probe. 

Generally,  the  following  measurements  were  taken  In  an 
experimental  run: 

1.  Relative  Irradlance  profile  across  the  laser  Image 
on  the  retina. 

2.  Horizontal  temperature  profile  in  the  axial  hot 
spot,  using  10  second  or  shorter  laser  pulses. 

3.  Temperature-time  history  for  a 10  second  or  shorter 
laser  pulse,  measured  at  the  center  of  the  image  and 
in  the  axial  position  of  highest  temperature  rise. 

4.  Axial  temperature  profile  through  the  pigment 
epithelium  and  choroid  following  1.6  millisecond 
to  10  second  laser  pulses. 

5.  Temperature  rise  at  the  center  of  the  image  versus 
corneal  power. 

A horizontal  scan  In  the  vitreous  humor  was  taken  for 
each  run  In  the  following  manner.  The  probe  was  Inserted  Into  the 
vitreous  where  It  was  clearly  visible  and  In  the  center  of  the 
retinal  Image.  The  animal  positioning  device  was  rotated,  which 
rotated  the  monkey  eye  and,  hence,  the  probe  tip  through  the 
retinal  Image.  The  retinal  Image  did  not  rotate  with  the  eye,  but 
remained  fixed  In  space  as  shown  schematically  In  Figure  6. 
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FIGURE  6.  Horizontal  Scan  of  Probe  Through  Retinal  Image 

The  probe  and  eye  were  rotated  In  steps  through  the 
retinal  Image.  At  each  step,  a 1.6  millisecond  pulse  of  light 
was  gated  into  the  eye  with  a shutter.  The  temperature  increase 
at  each  step  due  to  direct  light  absorption  In  the  thermocouple 
provided  a relative  measure  of  the  Irradlance  at  that  point  In  the 
plane  of  the  retina. 

As  described  previously,  a temperature  rise  due  to 
direct  light  absorption  was  readily  quantified  because  its  response 
time  was  much  faster  than  that  of  conducted  heat.  In  the  vitreous, 
the  separation  of  components  was  even  greater.  The  medium  sur- 
rounding the  probe  had  a very  low  absorption  Itself,  and  conduction 
of  heat  from  the  P.E.  required  roughly  a millisecond.  The  tempera- 
ture changes  due  to  direct  absorption  appeared  first  as  an  increase 
with  the  same  rise  time  as  the  corneal  power  pulse,  and  then  a 
step  decrease  of  the  same  magnitude  when  the  laser  was  cut  off. 

A calibration  relating  the  temperature  Increase  to  Irradlance 
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magnitude  was  not  performed;  the  values  obtained  were  relative 
ones. 

A micrometer  on  the  base  of  the  animal  positioning  device 


was  used  to  rotate  the  animal  and  provide  a calibrated  distance 
measurement  for  the  probe.  By  rotating  the  animal  about  the 
vertical  axis  of  the  eye,  the  laser  image  remained  fixed  in  space. 
Movement  of  the  probe  480um  In  the  plane  of  the  retina  corresponded 
to  approximately  a two-degree  rotation  of  the  eye.  The  image 
viewed  at  the  sclera  with  a fixed  microscope  did  not  appear  to 
move  In  space  when  the  animal  was  rotated. 

If  the  axis  of  rotation  were  in  front  of  the  cornea,  the 
laser  Image  would  follow  the  probe  and  artificially  broaden  the 
retinal  Image.  Likewise,  if  the  axis  were  behind  the  lens,  the  image 
would  be  artificially  narrowed.  It  should  be  pointed  out  that 
the  stability  of  the  laser  Image  could  not  be  determined  by  ob- 
serving It  through  the  fundus  camera.  Of  significance  was  the 
fact  that  the  probe  tip  had  to  move  horizontally  only  a millimeter 
or  less  to  scan  the  complete  Image.  Thus,  small  angular  rotations 
were  adequate. 

The  Irradlance  profile  and  retinal  Image  radii  for  a run 
were  determined  from  a graph  of  probe  potential  versus  horizontal 
distance,  the  potential  being  measured  at  the  completion  of  the 
shutter  opening.  Profiles  were  normalized  to  a value  of  1.0  at 
the  center  of  the  image,  and  image  radii  were  read  from  the  graph. 
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Subthreshold  temperature  measurements  are  made  with  the 
probe  centered  on  the  laser  image  and  in  the  axial  hot  spot.  Several 

experimental  runs  for  subthreshold  laser  exposures  are  made  for 
exposure  durations  from  .002  second  to  10  seconds  and  retinal 
Images  of  100  to  200pm.  Thresholds  are  then  determined  at  posi- 
tions 1/2  to  2 millimeters  on  either  side  of  the  probe  insertion 
site  for  various  times  of  exposure.  The  threshold  temperature 
is  then  extrapolated  from  the  subthreshold  measurements  to  the 
threshold  laser  power  value.  Laser  power  is  measured  for  each 
rur.  The  relative  temperature  rise  In  the  eye  is  measured  with 
respect  to  a 37°C  reference  junction  which  is  maintained  within 
± .05°C  by  a proportionally  controlled  oven.  The  system  bandwidth 
for  the  thermocouple  measurement  Is  D.C.  to  lOKHz.  Measurement 
error  of  the  system  is  estimated  to  be  ±.2°C. 

An  ophthalmoscoplcally  visible  threshold  lesion  when 

viewed  through  the  fundus  camera  is  a circular  gray  or  white 
region  with  radius  of  25  to  200pm,  depending  on  the  duration 
of  exposure.  Image  size,  and  time  elapsed  between  Irradiation  and 
observation.  The  threshold  temperature  actually  occurs,  not  at 
the  center  of  the  Image,  but  some  distance  from  the  center 
probably  near  the  radial  boundary  of  the  lesion.  Consequently, 
once  the  lesion  radius  Is  measured,  the  threshold  temperature  is 
the  maximum  temperature  In  the  experimental  temperature-time 
profile  at  the  distance  equal  to  the  lesion  radius.  For  a 
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comparison  of  experimental  results  with  the  damage  model,  the 
radius  Is  measured  with  the  fundus  camera  reticle. 

In  Figure  7,  the  laboratory  arrangement  of  the  ex- 
perimental instruments  is  illustrated.  The  Argon  laser  provides  a 

maximum  power  of  one  watt  at  the  488oX  wavelength  used  in  this  re- 
search. The  shutter  In  front  of  the  laser  is  used  to  select  the 
exposure  duration.  For  subthreshold  exposures  and  for  threshold 
determinations,  the  laser  power  is  controlled  to  within  4%  by 
neutral  density  filter  placed  either  between  the  laser  and  the 
shutter  or  in  front  of  the  shutter.  Optics  in  the  laser  path 
are  used  to  modify  the  laser  beam  by  providing  one  or  two  lens 
focus  of  the  beam  near  the  cornea  of  the  animal’s  eye.  A beam 
splitter  allows  the  simultaneous  viewing,  measurement  and  ex- 
posure of  the  fundus.  The  fundus  camera  is  used  to  observe  and 
photograph  the  retina  before  and  after  lesion  production.  The 
power  for  each  exposure  is  measured  by  the  radiometer  through 
the  beam  splitter.  Calibrations  of  the  amount  of  light  trans- 
mitted to  the  radiometer  and  reflected  into  the  animal's  eye  are 
performed  with  the  radiometer. 


FIGURE  7:  Laboratory  arrangement  of  Instruments. 


RESULTS 


Threshold  temperatures  were  measured  in  29  primate  eyes 
for  both  macular  and  paramacular  exposures.  The  center  temperature 
for  a threshold  lesion  was  determined  at  the  center  of  the  retinal 
image  and  the  axial  position  of  maximum  temperature  rise.  For 
some  of  the  runs,  lesion  diameters  were  visually  measured  and 
correlated  with  horizontal  temperature  scans  to  determine  the 
temperature  rise  at  the  lesion  diameter.  In  a small  number  of 
experiments,  axial  direct  absorption  profiles  were  taken  for  a 
measurement  of  tissue  absorption  coefficients. 

Lesion  threshold  temperatures  were  measured  for  Argon 
laser  exposures  from  .002second  to  lOseconds  in  duration.  Both 
macular  and  paramacular  insertion  sites  for  50  and  100pm  (average) 
1/2-power  radius  images  were  used  to  determine  threshold  lesion 
temperatures.  Temperatures  were  measured  at  several  subthreshold 
corneal  powers  and,  for  most  exposures,  threshold  temperatures 
were  extrapolated  from  a corneal  power  50-75%  of  threshold. 

Horizontal  and  axial  profiles  were  obtained  from  most 
eyes  for  the  comparison  of  the  experimental  temperature  profile 
with  a computer  solution  of  the  heat  conduction  equation  (40). 

The  general  experimental  protocol  Included  horizontal  beam  scans 
In  front  of  the  P.E.  to  determine  the  image  size  at  the  retina  of 


W.L 


,JWW 


55 


the  monkey  eye.  Horizontal  temperature  scans  were  accomplished  at 
several  times  and  corneal  powers  with  the  thermocouple  placed  In 
the  hottest  position  In  the  beam  profile  and  tissue.  Axial  scans 
from  the  neural  layers  back  to  the  sclera  were  done  at  the  con- 
clusion of  the  experiment  for  a particular  Insertion. 

The  “threshold"  temperature  rise  In  this  report  Is  the 
temperature  rise  of  the  P.E.  In  the  center  of  the  beam  that  Is 
associated  with  the  appearance  of  an  ophthalmoscoplcally  visible 
lesion  5 minutes  post  exposure.  The  lesion  radius  "threshold" 
temperature  rise  Is  the  estimated  temperature  at  the  lesion 
radius  obtained  from  measurements  of  lesion  size,  lesion  center 
temperature,  and  horizontal  temperature  scans.  All  exposures 
presented  were  obtained  with  an  Argon  (488oJi  or  all  lines)  laser. 
Irradiation  duration  was  controlled  by  a mechanical  shutter  with 
an  opening  time  (for  the  laser  beam)  of  approximately  .0002  second. 

In  this  section,  the  threshold  measurements,  the  pre- 
diction of  temperature  in  the  retinal  tissues,  and  a comparison 
of  lesion  measurements  with  a simple  rate  process  model  will  be 
discussed. 

Temperature  Rise  Versus  Corneal  Power 

For  each  experiment,  several  subthreshold  corneal  powers 
are  used  to  construct  a temperature  versus  corneal  power  graph  for 
extrapolation  to  the  threshold  corneal  power.  In  Figure  8 Is  shown 


FIGURE  8:  Temperature  rise  versus  corneal  power  for 
five  experiments. 


a plot  of  temperature  rise  versus  corneal  power  for  five  experi- 
ments. The  straight  lines  through  the  data  are  drawn  with  a slope 
of  one  to  Illustrate  the  linearity  of  temperature  with  corneal 
power.  The  profiles  for  these  measurements  are  all  approximately 
truncated  gausslan  profiles  with  half  power  radii  ranging  from 
35pm  to  125pm.  The  corneal  power  Is  adjusted  with  neutral  density 
filters  and  the  temperatures  are  recorded  at  the  end  of  1 and  10 
second  laser  pulses. 

The  differences  between  the  center  temperatures  for  a given 
corneal  power  are  due  to  the  differences  in  half  power  radii  for 
the  experimental  temperature  versus  corneal  power  measurements. 

For  example,  the  profile  for  the  experiment  denoted  by  monkey  1 
has  a half  power  radius  of  35pm,  the  retina  profile  for  monkey  6 
has  a half  power  radius  of  98pm  and  the  half  power  radius  for  the 
retina  image  for  monkey  2 is  125pm. 

Threshold  Temperature  Measurements 

Threshold  temperatures  associated  with  the  appearance  of 
an  ophthalmoscoplcally  visible  lesion  5 minutes  post  exposure  are 
measured  both  at  the  center  of  the  image  and  at  the  lesion  radius. 
The  center  threshold  temperature  Is  the  Increase  In  temperature 
relative  to  37°C  baseline  fundus  temperature.  The  lesion  radius 
temperature  measurement  Is  determined  from  subthreshold  corneal 
power  and  horizontal  (radial)  scans  of  the  retinal  temperature 
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profile.  The  temperature  at  the  lesion  radius  Is  extrapolated 
from  the  horizontal  scan  and  the  ophthalmoscopic  measurement  of 
the  gray  or  white  lesion  radius.  The  measurements  of  temperature 
Increase  are  accomplished  with  the  thermocouple  probe  described 
earlier,  while  the  thresholds  are  determined  from  a series  of 
exposures  near  the  probe  Insertion  site  for  varying  corneal 
power.  The  corneal  power  for  the  threshold  is  adjusted  by  neutral 
density  filters  from  a value  50  to  100*  above  threshold  to  20  to 
50%  below  threshold.  The  value  of  corneal  power  which  produces  a 
burn  which  Is  visible  5 minutes  post  exposure  is  the  threshold 
corneal  power.  Corneal  powers  above  these  values  always  produce 
a lesion  before  5 minutes  has  elapsed.  Power  levels  below  the 
threshold  corneal  power  do  not  produce  a lesion  within  5 minutes. 

A suirmary  of  the  threshold  temperature  rise  at  the  center 
of  the  lesion  for  exposure  durations  from  .002  second  to  10  seconds 
is  presented  in  Table  V.  The  values  in  this  table  are  extrapolated 
from  subthreshold  temperature  measurements  at  the  center  of  the 
image  and  at  the  axial  position  of  highest  temperature  rise.  All 
values  In  the  table  are  referenced  to  the  baseline  fundus  tempera- 
ture of  37°C.  Both  macular  exposures  and  paramacular  exposures 
(temporal  to  the  macula)  are  presented.  The  number  of  thresholds 
determined  at  each  exposure  duration  with  the  mean  temperature 
rise  and  standard  deviation  are  listed. 
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TABLE  V:  Summary  of  Threshold  Temperature  Measurements 
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Macular 
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Mean 
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1.7 

13 
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39.9 
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29.4 

4.1 
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32.0 

.05 

1 

57.6 
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7.4 
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3.8 
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39.2 
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43.0 
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At  10  seconds,  7 paramacular  and  13  macular  threshold 
measurements  yield  average  temperature  increases  of  24.8°C  and 
19.9°C  respectively;  the  standard  deviations  of  the  10  second 
threshold  temperature  increases  are  1.7°C  for  paramacular  expo- 
sures and  1.5°C  for  macular  exposures.  For  one  second  exposures, 
the  temperature  rise  required  to  produce  a visible  lesion  In- 
creases to  28.7°C  and  27.1°C  for  paramacular  and  macular  insertion 
sites.  At  0.1  second  paramacular  and  macular  temperature 
Increases  are  39.9°C  and  29.4°C  respectively. 

Generally,  the  standard  deviation  of  the  temperature 
Increases  Is  less  for  longer  curatlon  exposures  than  It  Is  for 
shorter  duration  exposures.  Even  for  longer  time  measurements 
for  which  the  number  of  measurements  Is  nearly  the  same  as  for 
shorter  times  of  exposure,  the  standard  deviation  Is  smaller. 


Expressed  as  percentages  of  the  peak  temperature  the  paramacular 
standard  deviations  for  10  second,  1 second,  and  .1  second  duration 
are  7%,  10%,  and  13%,  respectively;  macular  temperature  Increase 
standard  deviations  are  8%,  14%  and  14%  for  the  10  second,  1 
second,  and  .1  second  exposures. 

In  Figure  9 is  illustrated  the  average  temperature  rise 
versus  exposure  duration.  The  figure  is  a plot  of  the  data  in 
Table  V.  The  temperature  increase  for  the  appearance  of  a visible 
lesion  is  lower  for  macular  exposures  than  for  paramacular  ex- 
posures. This  figure  illustrates  a temperature-time  dependence 
for  lesion  production  which  may  be  different  for  macular  and  para- 
macular sites.  For  an  exposure  duration  of  10  seconds  the  macular 
temperature  increase  is  80%  of  the  increase  required  to  produce 
a paramacular  lesion.  For  1 second  and  .1  second  exposures,  the 
macular  increase  is  84%  and  74%  of  the  paramacular  temperature 
increase  respectively,  while  at  .02  second,  for  a small  number 
of  exposures,  the  macular  increase  is  91%  of  the  paramacular 
temperature  rise. 

For  a number  of  the  exposures,  lesion  radius  measurements 
are  accomplished  with  a reticle  in  the  fundus  camera  eyepiece. 

In  Table  VI,  a summary  of  the  measurements  of  lesions  radius  and 
the  extrapolated  temperature  from  horizontal  profile  measurements 
is  presented.  The  number  of  lesion  measurements,  the  exposure 
duration,  the  lesion  radius  normalized  to  the  half  power  radius 
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of  the  retinal  image,  and  the  temperature  at  the  lesion  radius 
normalized  to  the  center  temperature  are  tabulated  for  exposure 
duration  of  .01  second  to  10  seconds.  At  10  seconds,  thirteen 
lesion  measurements  yield  an  average  lesion  radius  of  .55  times 
the  half  power  radius.  The  temperature  at  the  edge  of  the  lesion 
is  .88  times  the  center  temperature  from  a 10  second  horizontal 
scan.  The  values  in  the  table  show  that  the  normalized  average 
temperature  at  the  lesion  radius  is  lower  for  short  duration  ex- 
posures than  for  long  exposure  times.  Also,  the  average  value 
of  normalized  lesion  radius  Increases  with  increasing  exposure 


duration. 


i 


TABLE  VI:  Lesion  Radius  Measurement  Normalized 

to  One-Half  Power  Radiance  and  Lesion 
Temperature  Normalized  to  Center  Tem- 
perature for  Macular  Exposures 


Number  of  Lesions  Time  of 
Measured  Exposure 


Lesion  Rad. 

I /2-Power  Rad. 
(average) 


■>.  at  Lesion  Rad. 
Center  Temp, 
(average) 
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Thermal  Model  Parameters 


The  source  term  In  the  heat  conduction  equation  Includes 
the  absorption  of  the  various  layers  In  the  ocular  fundus  and  the 
retinal  distribution  of  the  Input  light  energy.  The  source  term 
also  depends  upon  the  value  of  the  transmission  of  the  preretinal 
ocular  media  (TOM).  Tissue  thermal  parameters  conductivity  (k) 
and  volumetric  specific  heat  (pc)  also  affect  the  solution  for 
retinal  temperature  increase. 

In  Figure  10  Is  Illustrated  the  effect  of  changing  the 
transmission  of  the  preretinal  ocular  media  in  the  model.  Chang- 
ing the  transmission  of  the  preretinal  ocular  media  from  .79  to 
.51  In  Figure  10,  results  In  a reduction  in  the  temperature  In- 
crease for  all  times  in  the  10  second  temperature-time  history. 

The  transmission  parameter  of  .51  is  65%  of  the  parameter  .79. 

The  temperature  versus  time  for  the  two  values  of  transmission 
exhibit  the  same  decrease  of  35%  for  the  decrease  in  TOM  from  .79 
to  .51. 


In  Figure  11,  model  calculation  of  temperature  rise 
versus  corneal  power  Is  Illustrated  for  times  of  .01  second,  .1 
second  and  1 second.  The  temperature  Increases  for  a corneal 
power  of  10  milliwatts  are  half  the  Increases  for  20  milliwatts 
corneal  power.  The  temperature  Increase  for  all  times  of  expo- 
sure Is  linear  with  corneal  power.  An  Increase  In  corneal  power 
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» FIGURE  10:  Model  calculated  temperature- time  history  for  two 

values  of  transmission  of  the  pre-retinal  ocular 
Mdla  (TOM).  Values  for  TOM  of  .79  and  .51  are 
used  for  calculations  of  temperature  as  a function 
of  time  for  the  profile  and  absorpiton  parameters 
r • shown  on  the  figure. 
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FIGURE  11:  Calculated  temperature  as  a function  of  corneal  power 
for  times  of  .01  second,  .1  second  and  1 second  during 
a 1 second  exposure. 


Is  accompanied  by  the  same  percentage  increase  in  temperature 
rise. 

In  addition  to  the  transmission  of  the  preretinal  ocular 
media  (TOM)  and  the  corneal  power  input  to  the  model,  the  ab- 
sorption parameters  for  the  various  layers  of  the  ocular  fundus 
are  important  in  the  solutions  of  the  conduction  equation.  As 
mentioned  in  the  background  section,  several  sets  of  parameters 
have  been  calculated  from  measurements  of  absorption  of  light  in 
the  eye  and  its  component  parts.  In  Table  VII,  three  parameter 
sets  are  listed.  These  three  sets  are  used  for  comparisons  with 
the  experimental  temperature  measurements  in  the  next  section. 
Parameter  set  1 Is  taken  from  values  supplied  by  Brooks  Air  Force 
Base,  School  of  Aerospace  Medicine  (59).  The  transmission  of  the 
preretinal  ocular  media  Is  taken  to  be  .51  and  is  from  the  "direct" 
transmission  measurements  by  Boettner  (9).  The  absorption  coef- 
ficients for  the  pigment  epithelium  (APE)  and  the  choroid  (ACH) 
are  also  calculated  from  Boettner 's  measurements  which  do  not 
Include  differential  measurements  for  the  P.E.  and  choroid.  The 
coefficient  APE  represents  80%  absorption  in  the  P.E.  which  is 
probably  too  high  based  upon  measurements  by  Geeraets  and  Berry 
(7,  8). 

Parameter  set  2 Is  derived  from  measurements  by  Geeraets 
and  Berry  (7,  8).  The  transmission  of  the  preretinal  ocular  media 
represents  a trade-off  between  the  total  transmission  value  from 
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Geeraets  and  Berry  and  the  direct  transmission  value  from 
Boettner.  The  coefficient  APE  represents  a 53%  absorption  of  light 
energy  In  the  P.E.  The  choroid  absorbs  25%  of  the  light  reaching 
the  P.E.  This  implies  that  22%  of  the  light  entering  the  P.E.  is 
transmitted  through  the  choroid  to  the  sclera. 

Parameter  set  3 is  a combination  of  parameter  sets  1 and 
2 with  an  assumed  absorption  of  33%  of  the  light  energy  reaching 
the  P.E.  In  that  layer  and  58%  absorption  In  the  choroid.  Of 
the  light  reaching  the  P.E.,  9%  Is  transmitted  through  the  choroid 
Into  the  sc Vera. 

In  Figure  12,  the  absorption  of  light  versus  depth 
from  the  front  of  the  pigment  epithelium  for  the  three  parameters 
sets  In  Table  VII  is  Illustrated.  Parameter  set  1 represents  a 
high  absorption  of  light  in  the  P.E.  Parameter  set  3 gives  a much 
larger  absorption  of  light  In  the  choroid,  while  parameter  set  2 
represents  a lower  absorption  In  the  P.E.  than  set  1,  and  a lower 
absorption  of  light  In  the  choroid  than  set  3. 


TABLE  VII: 

Absorption  Parameters 

for  Figure  12 

SET 

TOM 

APE 

ACH 

1 

.51 

1622 

189 

2 

.636 

637.6 

76.28 

3 

.65 

395 

195 

ABSORPTION  PROFILES 


The  axial  distance 
pigment  epithelium. 


Model  temperature-time  histories  for  a 1 second  exposure 
for  the  three  parameter  sets  are  Illustrated  In  Figure  13.  The 
temperature- time  histories  were  generated  using  the  experimentally 
measured  retinal  Intensity  profile  from  MKY8-74  and  corneal  power 
of  26  milliwatts.  The  third  parameter  set  exhibits  a longer  time 
constant  than  the  first  set  since  the  10  millisecond  temperature 
Is  approximately  75*  of  the  temperature  rise  for  set  1,  while 
the  1.0  second  temperature  Increases  are  almost  equal. 

The  associated  temperature  decays  for  the  temperature 
rise  In  Figure  13  at  the  end  of  the  1 second  exposure  are  shown 
on  Figure  14.  The  decays  for  sets  1,2,  and  3 indicate  slightly 
different  time  constants  with  the  time  constant  for  set  3 being 
the  longest.  The  decay  time  (time  to  .368  peak)  for  set  3 is 
.03  second.  For  set  1 the  decay  time  is  .018  second  and  for  set 

2 It  Is  approximately  .021  second. 

Horizontal  scans  during  the  temperature  rise  for  the  three 
sets  of  model  absorption  parameters  are  Illustrated  in  Figure  15. 

The  temperature  field  at  .01  second  and  1.0  second  for  a standard 
one  second  exposure  are  Illustrated  for  all  parameter  sets.  The 
differences  in  time  constants  for  the  sets  may  be  seen  In  the 
comparison  of  the  .01  second  and  1 second  horizontal  scans.  The 
1 second  scan  Indicates  that  the  peak  temperatures  for  sets  1 and 

3 are  nearly  the  same  while  at  .01  second  the  peak  temperatures 
for  sets  1 and  3 are  nearly  the  same  and  the  peak  temperatures 
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FIGURE  14:  Model  calculated  temperature  decay  following  1 second  rise 
In  Figure  13  for  the  three  parameter  sets  In  Table  VII. 
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FIGURE  15:  Comparison  of  model  calculated  horizontal 
scans  for  three  parameter  sets  for  a 1 
second  exposure. 
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for  sets  2 and  3 are  closely  matched.  The  1 second  profile 
for  set  3 is  spread  more  radially  than  set  1. 

The  ratios  of  peak  temperatures  at  .01  second  to  peak 
temperatures  at  1 second  are  .53,  .49,  and  .44  for  sets  1,  2,  and 
3,  respectively.  The  time  constant,  therefore,  increases  with 
decreasing  absorption  in  the  pigment  epithelium.  This  reduces 
the  P.E.  source  term  and  reduces  short  time  temperatures.  At 
long  times,  the  temperature  is  a function  of  the  energy  deposited 
in  the  P.E.  and  choroid  and  the  tissue  conductivity.  The  total 
energy  deposition  Is  approximately  equal  for  parameter  sets  1 and 
3,  whereas  the  total  energy  deposition  for  set  2 is  less  than  sets 
1 and  3. 

The  axial  profiles  shown  in  Figure  16  illustrate  the  dif- 
ferences in  model  tissue  volume  heated  by  the  energy  deposited  in 
retinal  tissue  by  the  three  model  parameter  sets.  The  third  set 
heats  a larger  volume  of  tissue  to  a higher  temperature  for  the 
same  irradlance  profile  at  the  front  of  the  P.E.  The  axial  posi- 
tion of  highest  temperature  rise  for  sets  1 and  2 Is  at  the  center 
of  the  P.E.,  while  for  set  3 It  is  approximately  20ym  behind  the 
front  of  the  P.E. 

Once  appropriate  values  have  been  chosen  for  the  trans- 
mission, absorption  and  thermal  parameters  In  the  model,  one  more 
Important  parameter  remains.  The  radius  of  the  retinal  Image  Is 
extremely  Important  In  the  accurate  prediction  of  temperature. 


FIGURE  16:  Model  calculated  axial  scans  for  three  parameter 
sets  for  a 1 second  exposure. 


For  small  Images,  the  rise  time  of  the  temperature  Increase  near 
the  center  of  the  Image  is  much  faster  than  the  rise  time  for  larger 
Images.  The  Image  radius  at  the  retina  is  measured  during  a hori- 
zontal scan  with  the  thermocouple  in  the  neural  layers  of  the 
retina  at  a point  where  a large  direct  absorption  and  small  tissue 
rise  Is  present.  Alternative  methods  for  the  determination  of  image 
size  at  the  retina  are  from  the  geometric  optics  of  the  imaging 
system  or  from  the  point  spread  function  of  the  eye  convolved  with 
the  geometric  Image. 

The  effect  of  image  radius  on  the  rise  time  of  computed 
temperature  time  histories  is  shown  in  Figure  17.  For  example, 
a 60%  error  In  the  measurement  of  maximum  image  radius  (RIM)  when 
the  actual  value  Is  160pm,  would  produce  errors  of  56%,  18%  and  3% 
at  .01  second,  .1  second,  and  1.0  second  in  a 10  second  exposure 
calculation  for  the  normalized  temperature  Increase.  However, 
measurement  errors  of  10%  in  RIM  would  not  produce  excessive 
errors  In  the  computations  of  temperature  for  times  greater  than 
.01  second.  Five  different  model  parameter  sets  are  used  for 
calculation  of  the  temperature  rise  in  this  figure. 

For  the  model  set  with  parameter  set  3,  the  time  to  90% 
of  the  10  second  temperature  for  a 23pm  image  Is  22  milliseconds 
which  Is  nearly  the  same  as  for  the  model  with  set  1 for  a 40pm 
Image  radius.  The  model  parameter  set  3 (395,  195,  .65)  exhibits 
a longer  rise  time  for  a 150pm  Image  than  model  set  1 (1622,  189, 


Huttos/ 


for  varying  Image  radii.  Five  model  parameter  sets  are 
used  to  calculate  the  temperature  rise.  Temperature  time 
profiles  are  normalized  to  the  10  second  temperature  rise 
for  the  particular  beam  profile  and  model  parameter  set. 


77 


.51)  for  a 160pm  Image  radius.  For  model  set  1 the  time  to  90S  of 
the  10  second  temperature  versus  Image  radius  Is  22  milliseconds 
for  a 40pm  Image  and  5.7  seconds  for  a 2000pm  Image  radius.  The 
time  to  90%  of  the  10  second  temperature  for  parameter  set  3 is 
.29  seconds  for  a 150pm  radius  and  22  milliseconds  for  a 23pm  Image 
radius.  The  ratio  of  the  heated  volumes  for  set  3 and  the  150pm 
Image  to  the  23pm  Image  Is  42.5  while  the  rise  time  ratio  Is  12.6. 

Comparison  of  Experimental  and 
Model  Threshold  Temperatures 

For  most  of  the  experimental  thresholds  the  solution  to 
the  heat  conduction  equation  Is  computed  with  the  measured  retinal 
profile  and  the  threshold  corneal  power  for  the  absorption  and 
transmission  parameter  sets  given  In  Table  VII.  The  data  are 
grouped  according  to  lesion  location  (macular  or  paramacular). 

Within  each  group  the  data  are  divided  by  the  corneal  power  required 
to  produce  a visible  lesion  for  a given  Image  size  and  the  general 
shape  of  the  Image.  A larger  image  requires  a higher  corneal  power 
to  produce  a threshold  lesion  than  a smaller  Image.  If  the  half- 
power points  are  nearly  equal  the  tall  of  the  dlstrlt  .tlon  of 
retinal  power  Is  larger  for  one  Image  than  another,  the  Image  with 
large  tails  requires  a higher  power  to  produce  a visible  lesion. 

In  Table  VIII,  threshold  lesion  temperatures  at  the  radial 
center  of  the  Image  and  the  axial  position  of  highest  temperature 


...  it  \ ■ 


Increase  are  presented  for  a 10  second  exposure.  In  the  table, 
the  maximum  image  radius  (RIM),  the  half-power  radius,  experimental 
and  model  temperatures  for  macular  (M)  and  paramacular  (P)  inser- 
tion sites  and  the  threshold  corneal  power  are  presented  for  each 
experiment.  For  13  macular  Insertions,  the  mean  experimental 
temperature  rise  Is  19.9°C  with  a standard  deviation  of  1.5<>C. 

For  model  set  3 the  average  temperature  Increases  for  macular  and 
paramacular  exposures  are  20.4  ± 2.4°C  and  24.9  ± 5.3°C  respectively. 
No  average  values  for  corneal  power  are  given,  since  the  threshold 
corneal  power  depends  not  only  upon  the  image  size  and  distribution 
but  also  the  insertion  site. 

Threshold  data  for  one  second  exposures  are  shown  In 
Table  IX.  Experimentally  measured  temperatures  at  the  center  of 
the  visible  lesion  are  24.1°C  ± 3.4°C  for  macular  exposures  and 
28.7°C  ± 3.0°C  for  exposures  in  the  paramacular  region  temporal 
to  the  macula.  The  average  temperature  for  the  macular  exposures 
with  set  2 excludes  one  data  point  for  experiment  2/11/75  which 
Indicates  a much  higher  than  average  temperature  rise  of  60.7°C. 

Exposures  of  .1  second  duration  which  produce  threshold 
visible  lesions  and  the  associated  temperature  Increase  at  lesion 
center  are  listed  in  Table  X.  Fifteen  macular  exposures  and  8 
paramacular  exposures  yield  average  experimental  temperature 
increases  of  29.4°C  ± 3.9°C  and  40.0°C  ± 5.1°C.  The  total  corneal 
power  required  to  produce  a visible  lesion  5 minutes  post  exposure 
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Is  from  a minimum  of  28  milliwatts  for  a 152pm  maximum  radius 
(40pm  half-power)  radius  Image  to  105  milliwatts  for  a 480pm 
(98pm  half-power)  radius  Image,  (excluding  the  point  not  used 
In  temperature  average)  for  macular  exposures. 

Threshold  data  for  .01  second  exposures  are  presented  In 
Table  XI.  Experimental  measurements  at  3 macular  ar.d  4 paramacular 
Insertion  sites  yield  average  values  of  66.1°C  and  55.9°C  for  the 
threshold  lesion  center  temperature  increase.  For  the  paramacular 
exposures,  the  threshold  corneal  power  varies  from  90  milliwatts 
for  an  88pm  radius  Image  with  a half-power  radius  of  30pm  to  270 
milliwatts  for  a 248pm  radius  88pm  half-power  radius  image.  The 
macular  temperatures  probably  represent  some  anomalies  in  the  ex- 
periment since  the  corneal  power  for  a 72pm  image  does  not  differ 
significantly  from  the  corneal  power  for  a much  larger  285pm  image. 
For  these  two  images,  the  corneal  power  values  are  248  milliwatts 
and  214  milliwatts,  respectively. 

Miscellaneous  threshold  tenperatures  for  exposure  dura- 
tions from  .002  second  to  .05  second  are  present  in  Table  XII. 

The  data  Include  4 macular  exposures  of  .04  second  duration  with 
an  average  value  of  34°C  and  a standard  deviation  of  4.3°C.  For 
2 extramacular  Insertions  with  .04  second  duration,  the  average 
temperature  Increase  associated  with  the  appearance  of  a visible 
lesion  Is  57.8°C. 
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The  comparisons  of  model  and  experimental  threshold  tem- 
perature rise  versus  exposure  duration  are  sumnarlzed  In  Figures 
18  and  19.  In  Figure  18,  average  macular  experimental  threshold 
temperature  Increases  are  compared  with  the  solutions  of  the 
temperature  model  with  parameter  sets  1,  2,  and  3.  The  model  pre- 
dicts the  macular  tenperature  within  30%  at  all  times  of  exposure 
from  .04  second  to  10  seconds  for  parameters  set  1.  The  percentage 
errors  for  .04  second,  .1  second,  1 second  and  10  second  exposures 
using  set  1 In  the  temperature  model  are  +30%,  -.34%,  -9.5%,  and 
-3%,  respectively.  The  errors  for  exposure  durations  from  .04 
second  to  10  seconds  for  model  set  2 are  +14%  at  .04  second,  +1.4% 
at  .1  second,  +3.3%  at  1 second  and  -8%  at  10  seconds.  Errors  In 
the  prediction  of  threshold  temperature  using  the  model  with 
parameter  set  3 are  +2.7%  at  .04  second,  +9%  at  .1  second,  and  +1% 
for  1 second  and  +3%  at  10  seconds  exposure  duration. 

In  Figure  19,  extramacular  measured  and  predicted 
temperature  Increase  for  exposure  durations  from  .01  second  to 
10  seconds  are  displayed.  At  .01  second,  the  model  predicted 
temperature  Is  57%  higher  than  the  experimentally  determined  tem- 
perature Increase.  For  a .1  second  duration  Irradiation,  the  model 
predicts  a temperature  Increase  25%  higher  than  the  measured  rise. 
For  exposure  times  of  1.0  second  and  10  seconds,  the  model  predicted 
temperatures  with  parameter  set  1 are,  respectively,  4.9%  and  11.7% 
below  the  empirically  determined  temperature  Increases. 
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EXPOSURE  DURATION  IN  SECONOS 


FIGURE  19:  Comparison  of  experimental  temperature  Increases 
required  to  produce  visible  lesions  In  the  para- 
macuia  with  computed  values  from  the  conduction 
model  with  parameter  sets  1.  2 and  3. 
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Temperature-Time  Histories 

The  temperature  rise  required  to  produce  a visible 
lesion  is  a function  of  the  duration  of  the  irradiation,  the 
geometry  and  Intensity  within  the  image.  For  large  images,  the 
rise  time  of  the  temperature  at  the  center  of  the  Image  is  larger 
than  for  small  images  as  demonstrated  in  a previous  section. 

Also,  macular  and  paramacular  temperature  Increases  for  the 
appearance  of  a visible  lesion  are  clearly  different.  For  these 
reasons  and  for 'the  purpose  of  explanation  of  the  experimental 
factors  in  temperature  measurement,  this  section  will  deal  with 
typical  temperature- time  histories  for  horizontal  and  axial  scans 
In  macular  and  paramacular  Insertion  areas. 

In  Figure  20,  normalized  temperature-time  histories  for 
50um  and  100pm  half-power  macular  images  and  for  50ym  half-power 
paramacular  images  are  illustrated.  The  temperature  versus  time 
histories  represent  average  fractional  temperature  rise  (as  a frac- 
tion of  the  maximum  temperature  achieved)  versus  time  for  a 10 
second  laser  exposure  duration.  The  average  rise  time  to  .632  of 
the  10  second  temperature  for  a macular  50pm  half-power  radius  Is 
77.6  milliseconds.  The  100pm  half-power  macular  Images  require  an 
average  of  176  milliseconds  to  reach  .632  of  the  10  second  tempera- 
ture Increase.  A 50pm  half-power  radius  paramacular  temperature 
Increase  requires  a time  of  228  milliseconds  to  reach  .632  of  the 
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10  second  Increase. 

Temperature-time  histories  for  monkey  11-73  are  plotted 
In  Figures  21,  22,  and  23.  In  Figure  21,  an  experimental  tem- 
perature-time history  at  the  center  of  a horizontal  scan  Is  com- 
pared to  model  calculated  temperature  histories  for  the  parameter 
sets  1,2,  and  3.  At  .01  seconds,  the  model  (parameter  set  1) 
predicts  a temperature  Increase  of  6.1°C  while  the  measured  value 
Is  2.15°C.  Parameter  sets  2 and  3 cause  the  prediction  of  tem- 
perature at  .01  seconds  to  decrease  to  3.9°C  and  4.3°C  respectively. 
For  this  run,  the  measured  beam  profile  has  a maximum  radial  extent 
or  40ym  and  a half-power  radius  of  75ym.  Comparison  of  the  rise 
time  for  this  experiment  shows  that  the  rise  time  is  longer  than 
the  average  macular  rise  time  for  100ym  half-power  images.  The 
reason  for  this  could  be  the  large  tails  under  beam  profile. 

Figures  22  and  23  Illustrate  the  differences  in  rise  time 
for  two  different  positions  within  a horizontal  scan.  The  tem- 
perature-time history  In  Figure  21  Is  taken  from  Figure  22,  which 
Is  the  temperature  rise  versus  time  for  a 10  second  laser  irradia- 

i 

tlon.  The  rise  time  to  .632  of  the  10  second  temperature  In  Figure 
22  Is  approximately  .4  seconds.  The  oscillations  in  the  temperature 
time  history  are  caused  by  oscillations  of  ±3%  In  the  laser  output 
power.  The  decay  time  to  .368  of  the  10  second  temperature  for 
the  temperature  rise  at  the  center  of  the  horizontal  profile 
(Figure  22)  is  .42  second  which  is  nearly  the  same  as  the  rise  time 


FIGURE  21:  Experimental  and  model  calculated  temperature- time  histories 
for  monkey  11-73  for  a 10  second  exposure  duration.  Model 
calculations  of  temperature  Increase  with  three  parameter 
sets  are  plotted  with  measured  values. 
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during  the  laser  Irradiation.  The  direct  absorption  for  this  run 
Is  not  measurable  above  system  noise. 

In  Figure  23,  the  temperature- time  history  for  a position 
160*111  from  the  center  of  the  beam  profile  (85*jn  greater  than  the 
half-power  radius)  Is  shown.  The  profile  is  plotted  from  digitized 
data  on  the  Calcomp  plotter.  The  rise  time  at  this  position  off 
the  center  of  the  image  increases  to  1.6  seconds.  The  decay  time 
is  no  longer  nearly  equal  to  the  rise  time  for  temperature  rise  off 
center.  The  decay  time  for  this  position  has  increased  to  3 
seconds  due  to  conducted  heat  from  the  center  of  the  image. 

In  Figures  24  through  27,  the  effect  of  longer  paramacular 
rise  time  than  macular  rise  time  on  the  agreement  between  model  and 
experimental  temperature-time  histories  is  illustrated.  In  Figure 
24,  a .1  second  paramacular  temperature  rise  is  illustrated  for 
monkey  8-74.  Model  values  for  parameter  sets  1 and  2 are  plotted 
on  the  figure  for  comparison.  The  temperature  increase  for  the 
50pm  half-power  image  (experimental)  follows  a more  nearly  ramp- 
like increase  in  temperature  as  a function  of  time,  while  both 
models  exhibit  a much  faster  rise  time.  The  model  calculated 
tempera ture-time  history  could  be  approximated  more  closely  by  a 
step  increase  in  temperature  rise. 

The  variation  in  the  rise  time  of  temperature  for  positions 
±200ym  off  the  image  center  and  at  the  center  of  the  image  are 
illustrated  in  Figures  25,  26,  and  27.  In  Figures  25  and  27,  the 
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FIGURE  24:  Experiments 
time  hlstoi 
Is  linear. 


at  a position  200pm  to  the  observer's  left  of  bei 
In  the  axial  position  of  highest  temperature  rl s< 


center.  The  probe  Is  In  the  axial  position  of  hlghe: 
temperature  rise. 
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temperature  increase  is  nearly  linear  with  time  during  the  laser 
exposure.  The  maximum  temperatures  for  these  figures  are  3.0°C  and 
2.8°C,  respectively.  At  the  center  of  the  Image,  the  rise  of  tem- 
perature is  not  as  ramp-like  as  the  increase  200pm  off  center.  The 
maximum  Increase  is  9.3°C.  Near  the  center,  the  decay  of  the  tem- 
perature Is  more  rapid  relative  to  the  peak  temperature  than  at 
two  positions  off  center  due  to  heat  being  conducted  away  from 
the  center  and  toward  the  positions  off  the  image  center. 

In  Figure  28,  model  and  experimental  temperature  rises 
for  a 1 second  exposure  in  the  macula  of  monkey  4-74  are  illustrated. 
Model  predicted  temperature  increases  for  parameter  sets  1 and  3 are 
shown  in  the  figure  with  an  experimental  temperature  rise  at  the 
center  of  a 237pm  maximum  radius  image  with  a half-power  radius 
of  80pm.  The  model  and  experimental  values  match  closely  for  times 
during  the  exposure  from  .001  second  to  1 second.  The  rise  time 
(assuming  the  1 second  temperature  is  .87  times  the  10  second 
temperature  Increase)  is  approximately  200  milliseconds  which 
compares  favorably  with  the  rise  time  of  176  milliseconds  JFor  the 
"average"  100pm  half-power  Image  macula  temperature- time  history. 

In  the  measurement  of  temperature  Increase  and  the 
location  of  the  axial  position  of  maximum  temperature  rise,  the 
probe  Is  withdrawn  In  small  steps  from  the  vitreous  humor  Into  the 
retina  and  to  the  pigment  epithelium  (theoretically  the  position  of 
highest  temperature).  All  of  the  temperature- time  histories 
presented  In  this  section  were  taken  at  the  axial  position  of  highest 
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FIGURE  28:  Temperature  versus  time  for  a 1 second  exposure 
In  the  macula  of  monkey  4-74.  The  experimentally 
measured  temperature  rise  Is  plotted  for  comparison 
with  the  calculated  rise  for  model  sets  1 and  3. 
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temperature  rise.  Generally,  two  components  of  temperature  rise 
may  be  observed  on  the  thermocouple  output,  direct  absorption  and 
temperature  rise.  For  low  power  runs,  however,  the  direct  absor- 
ption component  may  be  absent  (lower  than  system  noise). 

In  Figures  29  through  34,  single  temperature- time  histories 
from  an  axial  scan  from  moneky  8-74  are  Illustrated.  The  probe  Is 
first  peaked  horizontally,  vertically  and  axially  In  the  position 
of  highest  temperature.  The  thermocouple  is  then  driven  axially 
with  the  hydraulic  drive  on  the  micromanipulator,  360pm  in  front 
of  the  P.E.  at  an  insertion  angle  of  10°.  Withdrawal  of  the  probe 
Is  accomplished  In  10pm  steps  until  the  P.E.  is  reached  and  then 
In  5pm  steps  for  the  next  120pm.  Ten  micron  steps  are  then  taken 
for  the  remainder  of  the  scan.  In  Figure  30  is  a computer  plot  of 
the  digitized  data  for  a position  280pm  in  front  of  the  P.E.  A 
direct  absorption  component  of  1.8°C  and  a temperature  rise  of  .7°C, 
which  may  be  a combination  of  tissue  rise  and  conduction  from  the  hot 
spot  down  the  probe.  Is  seen  on  the  record.  The  pulse  duration  for 
the  laser  Irradiation  Is  .01  second. 

A temperature  rise  and  direct  absorption  record  for  a 
probe  placement  80pm  In  front  of  the  P.E.  Is  shown  In  Figure  30. 

A direct  absorption  Increase  of  4.8°C  and  a temperature  rise  due  to 
conducted  heat  of  3.6°C  Is  present  on  the  record.  In  Figure  31,  a 
temperature- time  history  10pm  In  front  of  the  axial  hot  spot 
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Indicates  a direct  absorption  component  of  4.7°C  and  a conduction 
component  of  5.9°C.  For  the  10  millisecond  exposure.  In  the  hot 
spot,  the  conduction  temperature  rise  has  Increased  to  6.0°C  and 
the  direct  absorption  component  Is  4.6°C.  Two  other  runs  20pm 
and  60pm  behind  the  P.E.  are  shown  In  Figure  33  and  34.  At  the 
position  20pm  behind  the  P.E.,  the  direct  absorption  component  is 
3.5°C  and  the  temperature  Increase  due  to  conduction  Is  7.6°C.  At 
the  position  60pm  behind  the  P.E.,  the  direct  component  is  3.2°C 
and  the  temperature  Increase  Is  7.2°C. 

Horizontal  Temperature  Scans 

Horizontal  temperature  scans  are  performed  by  rotating 
the  animal  about  a center  of  rotation  near  the  nodal  point  of 
the  eye.  The  effect  of  this  rotation  is  to  move  the  image  past 
the  thermocouple  probe  and  to  thus  effect  a measurement  of 
temperature  rise  as  a function  of  horizontal  position  for  a given 
image  distribution.  Scans  obtained  in  this  manner  are  compared 
to  temperature  model  solutions  of  the  heat  conduction  equation. 

In  Figure  35,  a horizontal  (radial)  scan  for  monkey  11-73 
Is  Illustrated.  This  scan  Is  obtained  for  the  same  run  that  Is 
shown  In  Figure  21  In  the  section  on  temperature- time  histories. 

The  retinal  Image  distribution  for  this  macular  insertion  has  a 
maximum  Image  radius  of  400pm  and  a half-power  radius  of  75pm. 
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Horizontal  scan  for  monkey  11-74  In  the 
macula.  The  10  second  exposure  duration 
scan  Is  compared  to  the  model  solution  with 
parameter  set  1 at  various  times  during  the 
exposure. 


A large  amount  of  energy  is  therefore  present  in  the  tails  of  the 
distribution.  Measured  temperatures  at  times  of  .1  second,  1 
second,  and  10  seconds  during  a 10  second  duration  exposure  are 
compared  to  the  temperature  model  solutions  at  these  times  with 
the  absorption  parameters  of  set  1.  At  10  seconds,  at  the  center 
of  the  image,  the  model  predicts  a temperature  Increase  of  15.3°C 
while  the  measured  temperature  use  Is  16.2°C.  As  the  probe  is 
scanned  horizontally  from  the  center,  the  10  second  model  predic- 
tions are  lower  than  the  measured  temperatures  at  all  radial 
positions.  At  1 second  and  .1  second  the  model  predicts  higher 
increases  In  temperature  than  the  experimentally  measured  values. 
This  profile  represents  an  average  scan  for  large-tailed  images. 
The  major  effect  of  the  large  tails  Is  to  Increase  the  rise  time 
and  therefore  the  time  to  achieve  stea<Jy-state  temperatures. 

A .1  second  horizontal  scan  generated  by  computer  from 
digitized  data  Is  Illustrated  In  Figure  36.  The  Insertion  area 
Is  extramacular  temporal  and  superior  to  the  macula  approximately 
2mm.  The  temperature-time  history  corresponding  to  this  scan  is 
In  Figure  24.  Temperature  Increases  at  times  of  .001  second,  .01 
second,  .05  second  and  .1  second  during  a .1  second  exposure  are 
shown  in  Figure  37.  The  Image  for  this  run  has  a 152jjm  maximum 
radius  and  a half-power  radius  of  50ym.  The  slow  rise,  as 
evidenced  by  the  small  temperature  rise  at  .01  second  relative  to 
the  .1  second  temperature  increase,  Is  typical  of  paramacular 
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FIGURE  36:  Experimental  horizontal  scan  for  a .1  second  exposure  In 
the  paramacular  region  of  the  eye  of  monkey  8-74.  Tempera 
ture  increases  at  times  of  .01  second*  .05  second  and  .01 
second  are  plotted  versus  ladial  position. 


Insertions.  At  .01  second,  at  the  center  of  the  profile,  the 
temperature  Is  26%  of  the  .1  second  temperature.  The  model  with 
paramater  set  2,  predicts  a ratio  of  the  .01  second  to  the  .1 
second  temperature  of  84%. 

In  Figures  37  and  38,  temperature  Increase  versus  hori- 
zontal position  is  illustrated  for  times  during  and  following  a 
.01  second  exposure.  For  these  runs,  the  direct  absorption 
component  of  the  temperature  Increase  has  not  been  removed.  The 
Image  distribution  for  these  runs  is  the  same  as  for  the  previous 
figure.  This  run  Includes  a direct  absorption  component  at  the 
center  of  the  beam  of  4.8°C.  After  removal  of  the  direct  absorp- 
tion component,  the  decay  time  in  Figure  24  is  .062  second. 

A horizontal  scan  for  a paramacular  insertion  from 
monkey  6-74  is  shown  in  Figure  39.  Scans  for  times  of  .01  second 
and  .1  second  during  a .1  second  duration  exposure  are  compared 
to  model  calculations  using  parameter  sets  1 and  3.  The  Image  used 
in  this  run  has  a maximum  radius  of  144pm  and  a half-power  radius 
of  75pm.  At  .01  second  and  .1  second,  both  models  predict  higher 
temperature  Increases  than  the  experimentally  measured  Increase. 

In  most  paramacular  runs,  the  temperature  rise  at  large  horizontal 
displacements  is  higher  than  the  model  Indicates.  At  the  center, 
as  noted  earlier,  the  rise  time  of  temperature  for  the  paramacular 
insertions  is  much  longer  than  the  model  rise  time. 


In  Figure  40,  a horizontal  profile  for  monkey  4-74  Is 
displayed.  At  times  of  .01  second,  .1  second  and  1.0  second, 
model  set  1 and  3 are  plotted  for  comparison  with  experimental 
measurements.  The  Image  for  this  run  Is  237pm  In  maximum  radial 
extent  with  a half-power  radius  of  80pm.  This  profile  is  very 
nearly  the  average  macular  profile.  At  the  center  of  the  image 
and  for  all  horizontal  positions,  the  model  predicted  temperatures 
are  lower  than  the  experimentally  measured  temperatures  at  1 second. 
For  the  times  of  .1  second  and  .01  second,  the  temperature  rise 
predicted  by  the  model  with  parameter  sets  1 and  3 are  higher 
than  the  experimentally  measured  temperatures  at  all  radial 
displacements  from  the  center  of  the  Image. 

A final  feature  of  horizontal  temperatures  scans  is 
Illustrated  In  Figure  41,  with  the  aid  of  a .1  second  scan  for 
monkey  10-74.  The  input  profile  to  the  model  is  circularly 
symnetrlc,  so  unsymmetric  profiles  are  averaged  by  displacement 
from  the  center  of  the  image.  The  model  profile  In  the  figure 
reflects  the  symnetry  of  the  Input  Image  distribution.  The  scan 
from  the  experimental  data  exhibits  the  effect  of  an  unsymmetric 
profile  In  the  real  situation.  The  resulting  prediction  for 
the  "idealized"  image  does  not  match  the  experimentally  measured 
temperatures  on  either  side  of  the  center  since  It  does  not 
utilize  the  actual  unsymnetrlcal  beam  distribution  at  the  retina. 
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Axial  Temperature  Scans 


Axial  temperature  scans  indicate  the  temperature  as  a 
function  of  position  In  the  various  layers  and  fluids  of  the 
eye.  Axial  scans  start  In  the  vitreous,  progress  through  the 
neural  retina,  through  the  P.E.  and  then  into  the  choroid. 

In  the  layers  and  fluids  in  front  of  the  P.E.,  the  temperature 
rise  recorded  by  the  thermocouple  is  due  to  direct  absorption 
of  light  energy  In  the  metal  films  of  the  thermocouple  and  a 
conducted  component  which  is  a combination  of  heat  conduction 
down  the  probe  from  the  hot  spot  and  tissue  temperature  rise  due 
to  conducted  heat.  The  direct  absorption  component  is  used 
to  determine  absorption  profiles  in  the  tissue  and  the  temperature 
rise  component  Is  compared  to  the  model  predicted  temperature 
rise  for  each  scan. 

An  axial  scan  from  360pm  in  front  of  the  P.E.  to  300pm 
behind  the  P.E.  for  monkey  8-74  is  shown  in  Figure  42.  The 
scan  is  taken  In  10pm  steps  from  360pm  In  front  of  the  P.E. 
to  the  front  of  the  P.E.,  then  In  5pm  steps  from  0 to  120pm 
behind  the  P.E..  From  120pm  to  300pm  behind  the  P.E.,  the  scan 
Is  again  taken  In  10pm  axial  steps.  The  probe  track  angle  for 
this  experiment  Is  approximately  10°.  At  a position  360pm  In 
front  of  the  P.E.,  the  rise  In  temperature  Is  almost  all  direct 
absorption  as  evidenced  by  the  .003  second  and  .01  second 


exposure  duration.  Note  the  expanded  scale  rrom  a 
position  of  Oym  to  120pm  from  the  front  of  the  P.E. 


temperatures  being  nearly  equal.  The  initial  direct  absorption 
temperature  rise  is  1.9°C.  As  the  probe  is  withdrawn  toward  the 
P.E.,  the  direct  absorption  component  increase  due  to  the  angle 
of  insertion.  At  the  P.E.,  the  direct  absorption  rise  Is  4.9°C. 

If  the  direct  absorption  component  Is  corrected  by  the  angle  of 
the  insertion  and  the  beam  profile  values,  the  direct  absorption 
rise  at  the  -360ym  position  is  nearly  equal  to  the  direct  absorp- 
tion at  the  P.E.  The  direct  absorption  component  drops  off 
rapidly  as  the  thermocouple  is  withdrawn  behind  the  P.E.  The 
temperature  rise  reaches  a peak  approximately  20ym  behind  the  P.E. 

The  temperature  increase  as  a function  of  depth  behind 
the  P.E.  Is  plotted  in  Figure  43.  This  temperature  increase  is 
compared  to  solutions  of  the  temperature  model  using  parameter 
sets  1 and  3.  The  temperature  increase  is  derived  from  Figure  43 
by  subtracting  the  direct  absorption  component  at  each  axial 
position  from  the  total  temperature  rise.  The  exposure  time  for 
this  scan  Is  .01  second.  Both  the  model  predicted  scan  for  para- 
meter set  3 and  the  experimental  profile  peak  approximately  20vm 
behind  the  front  of  the  P.E.  The  model  using  parameter  set  1, 
however,  predicts  a peak  temperature  near  the  center  of  the  P.E. 

From  Figure  42,  the  direct  absorption  versus  depth  In 
the  axial  direction  relative  to  the  front  of  the  P.E.  Is  plotted 
In  Figure  44.  Two  other  profiles  from  monkey  4-74  and  monkey 
9-74  are  also  plotted  on  the  figure.  The  scan  from  monkey  8-74 


FIGURE  43:  Axial  temperature  scan  for  a .01  second  expsoure  for 
monkey  8-74.  Experimental  values  are  plotted  with 
calculated  temperature  Increases  for  model  sets  1 and  3 
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is  a paramacular  insertion  while  the  other  two  profiles  are  for 


macular  insertion  sites.  The  profile  from  monkey  8-74  shows  an 


absorption  of  32%  of  the  light  reaching  the  P.E.  in  that  layer 


and  almost  half  this  light  is  still  present  80ym  behind  the  front 


of  the  P.E.  The  two  macular  scans  indicate  that  nearly  50%  of 


the  light  energy  present  at  the  front  of  the  P.E.  is  absorbed  in 


the  lOym  thick  P.E.  and  only  30%  of  this  energy  is  present  80ym 


behind  the  front  of  the  P.E.  Calculated  absorption  coefficients 


from  the  scan  for  monkey  8-74  are  390.2  cm-1  for  the  P.E.  and 


102.6  cm-1  for  the  choroid.  For  monkey  4-74  and  monkey  9-74  scans, 


the  average  values  for  the  P.E.  and  choroid  absorption  coefficients 
are  648  cm"^  and  119.4  cnH  respectively. 


An  axial  scan  for  monkey  6-74  is  illustrated  in  Figures 


45  and  46.  In  Figure  45,  the  decay  of  temperature  for  a 91  milli- 


second exposure  is  shown  at  times  of  .901  second,  .111  second. 


.141  second,  and  .181  second.  In  Figure  46,  the  axial  scan  is 


plotted  again  for  comparison  with  model  sets  1 and  3.  Both  models 


predict  higher  temperatures  than  the  experimental  values  for  axial 


positions  near  the  P.E.  at  times  of  .01  second  and  .091  second. 

For  large  axial  displacements,  model  set  1 predicts  lower  tempera- 


ture Increase  than  the  experiment  indicates.  The  model  predicted 


temperature  Increase  versus  axial  position  for  parameter  set  3 


peaks  some  distance  behind  the  P.E.  while  the  calculated  temperature 
for  set  1 peaks  near  the  center  of  the  P.E. 


FIGURE  45:  Experimental  axial  scan  for  monl 
epithelium  Into  the  sclera.  Se' 
.091  second  exposure  are  plottei 
P.E. 


Comparison  with  Published  Thresholds 


The  data  obtained  In  this  study  are  grouped  by  half-power 
radii  for  comparison  with  the  thresholds  published  by  other  In- 
vestigators. Most  of  the  published  data  for  the  argon  wavelengths 
are  for  minimal  image  sizes.  Therefore,  the  larger  images  used  in 
this  study  should  require  higher  power  for  lesion  production. 

Threshold  corneal  powers  required  to  produce  a macular 
threshold  lesion  In  this  study  are  compared  to  data  from  Gibbons 
(15)  and  Clarke  (1)  in  Figure  47.  The  data  of  Gibbons  and  Clarke 
are  for  one  hour  threshold  criteria,  while  the  exposures  in  this 
study  are  for  threshold  lesion  appearance  5 minutes  post  exposure. 
Also,  Gibbon's  and  Clarke's  data  are  for  minimal  images  versus  the 
180pm  and  280pm  Image  radius  exposures  In  these  experiments.  For 
large  Images,  the  corneal  power  required  to  produce  a lesion  Is 
higher,  as  would  be  expected,  than  that  required  to  produce  a 
lesion  for  minimal  Image  sizes.  The  shorter  time  criterion  for 
these  experiments  (5  minutes)  should  also  contribute  to  higher 
corneal  power  for  lesion  appearance  than  the  one  hour  criterion. 

In  Figure  48,  the  threshold  corneal  powers  required  for 
lesions  for  extramacular  exposures  are  compared  to  the  values 
from  Vassllladls  (6)  and  Clarke  (1).  Both  15pm  and  100pm 
Images  from  Vassllladls  data  are  shown  In  the  figure  with  data 
for  180pm  and  280pm  Image  radii. 
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For  comparison,  the  center  temperature  that  would  be 
produced  by  the  threshold  powers  reported  by  Gibbons  (15)  and 
Vassllladls  (6)  Is  calculated.  First,  It  Is  necessary  to  assume 
a retinal  profile.  Vassllladls  reports  an  Image  radius  of  approxl- 
mately  15pm  and  a normally  distributed  laser  beam.  Assuming  the 
15im  represents  the  1/e2  point  on  the  distribution,  the  computed 
standard  deviation  of  the  retinal  Image  Is  7.5ym  and  the  tall  of 
the  distribution  extends  to  22.8pm.  This  profile  Is  also  applied 
to  the  minimal  Image  data  of  Gibbons.  The  computed  temperatures 
along  with  the  measured  average  temperature  threshol ds  are  shown 
In  Figure  49.  The  most  striking  feature  of  the  curve  Is  the 
.0135  second  minimal  Image  data  of  Vassllladls.  The  point  Is 
unexpectedly  low  when  compared  to  the  other  curves. 

The  100pm  data  of  Vassllladls  matches  the  measured 
temperatures  with  the  exception  of  the  .0135  second  point.  As 
might  be  expected,  the  minimal  Image  data  for  both  models  Indicate 
a higher  center  temperature  than  for  the  larger  Images  of 
Vasslliadls  and  this  research. 

A further  comparison  of  the  average  threshold  data  with 
the  data  of  Clarke  (1)  Is  shown  In  Figure  50.  The  average  threshold 
Irradlance  (w/cm2)  for  this  study  for  75pm  (average)  image  half- 
power radius  Is  plotted  for  comparison  with  the  data  of  Ham  et  al. 
(5)  for  10pm  and  100pm  disk  Images.  The  average  threshold 
Irradlance  is  calculated  from  the  model  with  a transmission  of  the 


Comparison  of  measured  temperatures  with  calculated  tempera 
tures  for  the  threshold  corneal  powers  from  the  data  of 
Vassllladls  (6)  and  Gibbons  (15). 


— Data  from  Ham  (5) 
x Average  Threshold 
Values  from  model 


TIME-SECONDS 


FIGURE  50:  Retinal  power  density  versus  time  of  exposure 
for  Image  radii  from  10um  to  800ym  from  Ham  (5) 
plotted  with  power  density  for  average  thres- 
holds for  .1  second*  1 second  and  10  second 
exposures. 
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preretinal  ocular  media  of  .636.  For  Ham's  study,  the  trans- 
mission of  the  preretinal  ocular  media  was  assumed  to  be  95%. 
The  data  from  Ham  et  al.  is  for  a white  light  xenon  arc  lamp 
with  a filtered  flat  spectrum  in  the  visible  region. 

Rate  Process  Model 


The  prediction  of  temperature  rise  in  the  retina  is  not 
sufficient  for  the  prediction  of  damage  due  to  high  intensity 
light.  Some  "effect"  function  must  be  used  in  conjunction  with 
the  temperature  model  to  indicate  injury  at  the  damage  temperature. 
The  rate  process  model  proposed  by  Henriques  (51)  was  intended 
for  the  purpose  of  prediction  of  porcine  epidermal  injury  as  a 
result  of  thermal  insult.  The  form  of  the  model  is  discussed 
In  an  earlier  section  of  the  background.  The  temperature  measure- 
ments in  this  study  that  result  in  damage  are  compared  to  those 
predicted  by  the  model  proposed  by  Henriques.  Rate  model  coef- 
ficients are  also  calculated  for  exposure  durations  from  .1 
second  to  10  seconds  for  the  temperature  history  at  the  lesion 
radius  In  this  research. 

Paramacular  exposure  center  temperatures  which  produced 
an  ophthalmoscoplcally  visible  lesion  are  plotted  In  Figure  51 
with  the  lesion  radius  temperatures  and  the  temperature  of  damage 
predicted  by  Henriques  rate  process  model  for  fl  = 1.  The 
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coefficients  in  the  rate  process  model  are  A * 3.1  x 1098  sec”1 
and  E = 150,000  cal /mole.  At  the  lesion  radius,  the  rate  process 
damage  Integral  in  conjunction  with  the  thermal  model  using  para- 
meter set  3 predicts  temperatures  of  29.5°C  for  a 10  millisecond 
exposure,  26.5°C  for  a .1  second  exposure  and  21.8°C  and  18°C  for 
1 second  and  10  second  exposures  respectively.  The  model  values 
are  computed  for  a gaussian  Image  profile  with  half-power  radius 

t.  I 

of  100pm.  For  the  .01  second,  .1  second,  1 second  and  10  second 
exposure  durations,  the  experimental  lesion  radius  temperature 
increases  are  39.4°C,  29.9°C,  23.5°C  and  21 .8°C, respectively. 

The  higher  paramacular  lesion  radius  temperatures  cause  the  rate 
process  model  to  predict  a larger  radius  than  the  radius  observed 
in  the  experiment. 

The  macular  temperature  increases  displayed  in  Figure  52 
with  the  rate  process  predicted  temperature  show  a different  trend 
than  the  paramacular  temperatures  in  the  previous  figure.  For 
the  .1  second,  1 second  and  10  second  exposure  durations,  the 
experimental  lesion  radius  temperature  Increases  are  83%,  91% 
and  97%  of  the  damage  temperature  predicted  by  the  rate  process 
model  In  conjunction  with  the  thermal  model.  For  the  .01  second 
exposure  duration,  the  experimental  temperature  Is  151%  of  the 
predicted  damage  temperature  from  the  rate  process  model. 

A comparison  of  the  coefficients  of  the  model  by  Henri ques 
(51)  for  paramacular  and  macular  temperatures  with  the  experimental 
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temperature  Increases  may  be  made  with  the  aid  of  Figures  53 
and  54.  The  rate  process  model  coefficients  are  derived  from 
longer  exposure  durations  than  those  In  the  present  study  and 
are  obtained  by  assuming  a step  Increase  In  temperature.  The 
model  coefficients  are  derived  from  a plot  of  in  t versus  1/T 
where  t Is  the  exposure  duration  and  1/T  Is  the  reciprocal  of 
temperature  In  degrees  Kelvin.  The  slope  of  this  plot  Is  E/R 
while  the  Intercept  Is  -AnA  In  the  equation: 

fl  = Ae  "E^T  t 


for  an  assumed  value  of  8 = 1. 


In  Figure  53,  graphs  of  An  t versus  the  reciprocal  of 
the  temperatures  at  the  center  of  the  lesion  for  both  paramacular 
and  macular  exposures  and  for  Henrlques  model  assuming  a step 
temperature  increase  is  illustrated.  The  paramacular  and  macular 
temperatures  at  the  lesion  center  are  both  greater  than  the  tempera- 
ture of  damage  (temperature  at  which  ft  ■ 1)  for  all  exposure 
durations.  At  the  shorter  durations,  however,  the  departure 
from  a straight  line  plot  for  An  t versus  1/T  Is  readily  apparent. 
The  temperature  rise  at  the  center  Is  not  a good  Indicator,  since 
the  area  of  damage  extends  radially  some  distance  from  the  center. 

In  Figure  54,  the  paramacular  and  macular  temperature 
Increases  at  the  center  are  multiplied  by  the  average  factor  for 
temperature  at  the  lesion  radius  (see  Table  VI)  and  plotted  on 
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(EXPOSURE  DURATION) 


PARAMACULAR 

MACULAR 

HENRIOUES 


.6  2.7  2.8  2.9  3.0  3.1 

l/T  IN  #K“I  x I0"3 

FIGURE  53:  Plot  of  An  (exposure  duration)  versus  the 

reciprocal  of  temperature  rise  at  the  center 
of  nacular  and  paramacular  lesions.  Henrlques 
model  predicted  temperatures  are  plotted  for 
comparison  (51). 
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the  in  t versus  1/T  scale  for  comparison  with  the  model  by  Henri ques. 
From  .1  second  to  10  seconds  ( n t s -2.3  and  2.3  respectively)  the 
data  for  macular  and  paramacular  lesions  may  be  nearly  approximated 
by  straight  lines.  At  .01  seconds,  the  data  again  depart  signi- 
ficantly from  the  straight  line.  Coefficients  may  be  derived 
from  the  two  straight  line  plots.  The  values  of  in  t and  1/T 
used  to  derive  the  rate  process  coefficients  are  given  In  Table 
XIII.  The  calculated  coefficients  in  the  rate  process  model  for 
paramacular  insertions  are  E/R  * 71,956  cal/mole  and  A s 2.55  x 10^4 
sec-1  versus  75,000  cal/mole  and  3.1  x 10®®  sec  1 for  Henriques 


model.  For  macular  exposures,  the  calculated  coefficients  are 
E/R  * 112,300  cal/mole  and  A * 8.4  x 10147  sec-1. 


TABLE  XIII:  Values  Used  to  Derive 

Rate  Process  Coefficients 


1/T  (Paramacular) 
3.014  x 10"3 


1/T  (Macular) 


-2.3026 


2.95  x 10 


3.012  x 10 


Probe  Response  Measurements 

The  rise  time  for  the  direct  absorption  of  the  thermo- 
couple in  different  media  is  measured  with  an  acousto-optic  modulator 
In  the  laser  beam.  The  rise  time  of  the  laser  light  pulse  Is  less 
than  5p  seconds.  Rise  times  for  the  thermocouple  In  air,  water 
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In  (EXPOSURE  DURATION) 


FI6URE  54:  Plot  of  in  (exposure  duration)  versus  the 
reciprocal  of  lesion  radius  temperature 
compared  to  Henriques  model  prediction  of 
damage  temperature  (51). 


and  vitreous  media  are  measured  for  the  step-like  Increase  In  light 
input.  The  laser  light  must  be  focused  to  a small  Image  In  order 
to  achieve  temperature  Increases  which  can  be  measured  above  the 
noise  of  the  system.  Rise  times  are  measured  for  10  to  20°C  rises 
above  the  temperature  of  the  medium  In  which  the  thermocouple  Is 
Imbedded.  For  a 19ym  diameter  probe,  the  rise  time  from  10%  to 
90%  of  the  final  temperature  Is  .4  milliseconds  In  a water  medium. 
This  value  compares  to  a computed  rise  time  using  equation  34  from 
Reed  (49)  for  the  step  response  of  the  probe  In  a heated  medium 
of  .3  milliseconds.  In  vitreous  humor,  the  probe  rise  time  Is  .7 
milliseconds  which  compares  to  Reed's  computed  rise  time  of  .8 
milliseconds  In  a biomedium.  In  air,  the  measured  10-90%  response 
time  of  the  probe  Is  12  milliseconds.  This  value  Is  much  lower 
than  the  computed  value  of  30  milliseconds  from  Reed's  model  (49). 

In  a small  number  of  experiments,  the  amplitude  versus 
frequency  plot  for  a few  thermocouples  Is  run  in  air  and  water 
media.  The  frequency  response  is  run  with  the  acousto-optic 
modulator  modulated  with  a sinusoidal  voltage.  Focusing  Is  again 
used  to  Increase  the  power  density  at  the  thermocouple  junction. 

The  amplitude  of  the  thermocouple  response  versus  the  frequency 
of  the  Incoming  light  Is  plotted  for  both  the  air  and  water  media 
for  the  same  power  density  at  the  junction.  The  frequency  response 
of  the  thermocouple  In  water  Is  10-20  times  the  frequency  response 
(3d8  frequency)  In  air.  The  amplitude  of  the  thermocouple  response  Is 
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larger  for  air  than  for  water.  If  the  frequency  response  (or 
rise  time)  in  water  10  times  the  response  in  air  then  the  amplitude 
of  the  temperature  in  air  is  10  times  the  amplitude  in  water. 

In  Table  XIV,  frequency  response  measurements  for  thermo- 
couple  probes  in  air  and  water  are  tabulated  for  three  probe 
diameters.  The  rise  time  indicated  is  computed  from  the  3dB  fre- 
quency (the  frequency  of  modulation  of  the  light-input  for  which 
the  voltage  response  is  .707  times  the  input  voltage).  The  rise 
time  Is  related  to  the  3dB  frequency  (f3dB)  by  * .SS/f^g- 

The  frequency  response  is  greater  for  smaller  (19um)  probes  than 
for  larger  probes.  The  computed  rise  times  from  the  frequency 
response  tests  are  higher  for  air  than  for  water.  For  a 19pm 
diameter  probe,  the  average  rise  time  is  8.6  milliseconds  In  air 
and  .35  milliseconds  in  water.  For  the  24pm  and  30pm  diameter 
probes,  the  average  rise  times  in  air  are  15.7  milliseconds  and 
22.4  milliseconds,  respectively.  With  the  24pm  and  30pm  probes  In 
water,  the  average  rise  times  are  .63  milliseconds  and  1.4  milli- 
seconds, respectively.  The  ratios  of  the  average  probe  response 
time  In  air  to  the  average  response  time  In  water  for  the  19pm, 
24pm,  and  30pm  probes  are  24.6,  24.9,  and  15.9  respectively. 


143 


TABLE  XIV:  Frequency  Response  and  Computed 

Rise  Time  for  Three  Probe  Diameters  In 
Air  and  Water  Media. 


Probe 

Diameter 

(urn) 

3dB  Frequency 
Air  Water 

(Hz)  (Hz) 

Rise  Time 

Air  Water 

(Seconds)  (Seconds) 

19 

36.8 

1000 

.0095 

.0035 

19 

35.0 

875 

.010 

.0004 

19 

55 

1170 

.0064 

.0003 

24 

25.8 

556 

.0136 

.00063 

24 

17.8 

520 

.0197 

.00067 

24 

25.5 

600 

.0137 

.00058 

30 

15.0 

243 

.0233 

.00144 

30 

16.4 

220 

.0214 

.0016 

30 

— 

292 

— 

.0012 
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DISCUSSION 


Threshold  temperatures  have  been  measured  with  thin  film 
copper-nickel  thermocouples  for  the  appearance  of  a visible  lesion 
In  the  primate  eye  following  exposure  to  laser  Irradiation. 
Temperature-time  histories  and  profiles  of  the  temperature  Increase 
In  the  retina  have  been  presented  In  the  previous  section.  The 
results  have  been  compared  to  a digital  computer  solution  of  the 
heat  conduction  equation  and  a simple  rate  process  model  based  on 
first  order  reaction  rates. 

Temperature  Measurement 

The  measurement  of  temperature  In  a living  tissue  with 
an  Invasive  technique  presents  several  problems  which  affect  the 
accuracy  of  the  measurement.  The  probe  Is  not  truly  authothermlc 
with  the  tissue  and  is  of  finite  size.  The  differences  In  thermal 
and  absorptive  properties  of  the  probe  and  tissue  cause  errors  In 
the  probe  measurement.  Ideally,  the  probe  should  have  Identical 
properties  with  the  tissue  and  be  In  Intimate  contact  with  the 
system  under  measurement.  Since  a complete  match  of  the  probe  and 
tissue  properties  Is  not  possible,  care  must  be  exercised  In  order 
to  minimize  the  effect  of  probe  and  tissue  property  differences  on 
the  measurement  of  tissue  temperature. 

In  this  research,  probe  diameters  from  10pm  to  20pm 
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were  employed  to  minimize  the  effect  of  the  size  of  the  probe  on 
the  measurement  of  tissue  temperature.  The  finite  size  of  the 
probe  causes  a shadowing  effect  for  the  light  entering  the  tissue 
since  the  probe  absorbs  or  reflects  all  the  light  which  impinges 
upon  It.  To  minimize  the  effect  of  shadowing,  or  removal  of  a 
portion  of  the  source  term,  the  Image  diameter  must  be  much  larger 
than  the  probe  diameter.  If  the  image  is  10  times  the  size  of  the 
probe,  only  1 or  2%  of  the  source  term  will  be  lost. 
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Another  problem  associated  with  probe  size  is  the  averag- 
ing of  temperature  across  the  junction.  Horizontal  or  radial 
averaging  may  be  reduced  by  using  larger  images  since  large  images 
tend  to  have  smaller  gradients  in  temperature  with  radial  position. 
The  ocular  fundus  contains  two  major  absorbing  layers,  the  P.E. 
and  choroid.  THE  P.E.  is  only  lOym  thick,  and  is  highly  absorbing, 
so  small  insertion  angles  of  10-30°  cause  some  averaging  to  occur. 

A simple  series  of  averages  from  the  temperature  model  with 
parameter  set  1,  Indicate  that  a 30°  angle  causes  the  errors 
Indicated  in  Table  XIV.  The  model  averages  are  taken  for  half- 
power radii  greater  than  75ym.  No  horizontal  average  Is  taken 
since  the  model  profile  has  no  variation  In  the  center  of  the  beam 
for  the  lOym  or  20ym  probe  radius.  The  average  is  a simple  axial 
average  over  the  depth  represented  by  D sin  6 where  D Is  the  probe 
diameter  and  6 Is  the  probe  Insertion  angle.  For  the  20ym 
diameter  probe,  the  average  Is  taken  over  a lOym  depth  for  a 30° 
Insertion  angle.  The  temperature  variation  with  depth  for  the 


average  Is  assumed  to  be  linear  over  the  small  depth  with 

(^hot  + ^cool)  the  average  value  of  tissue  temperature.  The 

2 

data  In  Table  XV  are  computed  from  the  model  values  of  monkey 
9-73,  11-73,  1-74,  4-74,  9-74,  10-74,  and  11-74.  The  values  for 
RIM  and  half-power  radii  are  given  in  Table  X.  The  errors  listed 
In  Table  XV  are  maximum  errors  for  the  probe  positioned  near  the 
center  of  the  P.E.  (axial  position  of  highest  temperature  for  model 


set  1). 


TABLE  XV:  Average  Error  Due  to  Axial 

Averaging  of  Temperature  for  Two 
Probe  Diameters 

% Error  for  a 20pm  % Error  for  a 40pm 


Diameter  Probe 
0 " s,  n) 

26,  3,  7 

8,  1.2,  7 

4,  .9,  7 

3,  .7,  6 


Diameter  Probe 
(1  - s,  n) 

46,  4,  7 

17,  2,  7 

7,  1.8,  7 

6,  1.6,  6 


Averaging  over  the  axial  length  of  the  probe  junction  causes  errors 
for  a 20pm  probe  of  8%  at  .01  second  and  3%  at  1 second.  For  a 
40pm  probe,  the  error  due  to  averaging  Is  estimated  to  be  17% 
at  .01  second  and  6%  at  1 second. 

The  averages  above  represent  a worst  case  average  with 
parameter  set  1,  since  the  absorption  parameter  for  the  P.E.  Implies 
that  80%  of  the  light  energy  reaching  the  layer  Is  absorbed  In  the 
P.E.  Realistically,  for  parameter  sets  2 and  3,  the  amount  of 
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averaging  Is  greatly  reduced.  In  fact,  for  the  20pm  or  smaller 
probes  used  In  this  study,  and  parameter  set  2,  the  average  tempera- 
tures across  the  probe  Is  within  5*  at  .01  second  and  for  set  3, 
the  temperature  Is  within  3*  at  .01  second.  The  typical  probe 
Insertion  Is  nearer  15°  than  the  30°  used  to  calculate  the  effect  of 
averaging.  Therefore,  probe  averaging  does  not  present  a significant 
problem  and  Is  probably  less  than  3%  for  most  exposure  durations 
In  the  present  study. 

Another  effect  of  finite  probe  size  is  the  axial  resist- 
ance to  heat  conduction  down  the  probe.  Larger  probes  have  a 
larger  effective  cross-sectional  area  and  thus  a lower  resistance 
to  heat  conduction.  The  resistance  to  heat  conduction  out  the 
leads  is  Inversely  proportional  to  the  conductivity  of  the  lead 
material  and  the  cross-sectional  area  of  the  lead  and  directly 
proportional  to  the  effective  length  of  the  lead.  In  the  present 
thermocouples,  the  small  cross-sectional  area  of  the  high  conduct- 
ivity metal  leads  and  the  low  conductivity  of  the  quartz  substrate 
tend  to  reduce  the  lead  resistance  effect.  The  effect  of  the  lead 
resistance  Is  to  cause  the  probe  to  measure  a lower  temperature 
than  actual  tissue  temperature.  The  magnitude  of  this  error  Is 
difficult  to  determine  because  It  also  depends  upon  the  resistance 
of  the  probe-tissue  Interface.  This  resistance  Is  due  to  a finite 
heat  transfer  coefficient  between  the  thermocouple  probe  and  the 
medium  In  which  the  temperature  Is  to  be  measured.  For  a higher 
junction  to  tissue  resistance  and  a lower  lead  conduction  reslst- 
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ance,  the  probe  measures  a temperature  lower  than  for  the  reverse 
situation. 

An  estimate  of  the  effect  of  the  probe  lead  resistance 
upon  the  measurement  of  steady-state  temperature  may  be  made  with 
the  aid  of  a simple  resistance-capacitance  analog  of  the  probe- 
tissue  system.  In  Figure  55  below,  the  analog  to  the  measure- 
ment of  step  temperature  Increase  is  illustrated. 


FIGURE  55:  Resistance  and  Capacitance  Analog  for  the 

Thermocouple  Probe 

For  the  measurements  by  Cain  and  Welch  (32),  and  the  model  by 
Reed  (49),  the  probe  is  Inserted  at  t * 0 into  a medium  at  constant 
temperature.  The  resistance  RJT  represents  the  junction  (or  probe) 
to  tissue  thermal  resistance,  while  represents  the  thermal 
resistance  of  the  probe  leads.  At  steady-state,  the  temperature 
of  the  junction  Is  given  by 

T1n\ 
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The  magnitude  of  steady  state  error  In  the  probe  measurement 
depends  upon  the  relative  sizes  of  Rj_  and  Rjj.  The  rise  time  of 
the  probe  In  this  situation  depends  upon  the  value  of  the  effective 
thermal  capacitance  of  the  probe  (C)  which  Is  given  by  pVcp  where 


pCp  Is  the  volumetric  specific  heat  of  the  probe  and  V Is  the  volume 
of  the  probe  heated.  The  thermal  resistances  In  parallel  form  an 
effective  charging  resistance  for  the  probe  thermal  capacitance. 

The  rise  time  is  given  by: 
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For  large  values  of  RL  the  rise  time  reduces  to  RjyC  which  is  an 
Idealized  rise  time  and  Is  the  case  for  Reed's  analysis.  If  RJT 
is  larger  than  the  ideal  thermal  resistance,  the  rise  time  increases. 
In  the  experimental  measurements,  this  resistance  appears  to  be 
larger  than  ideal  since  the  rise  times  measured  by  Cain  and  Welch 
(62)  in  water  are  longer  than  the  computed  rise  times. 

If  the  lead  resistance  (R|^)  is  less  than  about  10  times 
the  junction  to  tissue  resistance  (Rjj).  a significant  temperature 
drop  between  the  tissue  and  the  thermocouple  exists.  The  heat  flux 
out  the  thermocouple  leads  is  given  by: 


where 

AT.  * the  temperature  drop  from  the  tissue  to  the 
environment  (37°C). 

qL  ■ the  heat  flux  out  the  leads  cal/sec 
■ the  lead  resistance  ■ p- 
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the  effective  lead  length 
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k » the  thermal  conductivity 
A = the  cross-sectional  area  of  the  leads. 


The  heat  flux  from  the  tissue  to  the  junction  at  steady  state  is 
equal  to  the  heat  flux  out  the  thermocouple  leads.  For  a hori- 
zontal cylinder  (64),  the  heat  transfer  coefficient  from  a fluid 
medium  (such  as  water)  for  natural  convection  is  given  by: 

k V4 

hc  = .53  J (GrPr) 

where 


hc  ■ the  heat  transfer  coefficient  from  the  fluid 
to  the  probe. 

k = the  thermal  conductivity  of  the  fluid 
0 = the  diameter  of  the  cylinder 

Gr  = the  Grashof  number  * f (ATjt) 


Pr  » the  Prandtl  number. 

The  heat  flux  from  the  medium  to  the  probe  is  given  by: 

Ojt  * - hcta-TL)  atjt 

rjt 

5/4 

* ^(ATjt  )• 

The  properties  of  the  medium  and  the  heat  fluxes  are  functions 
of  temperature  and  temperature  driving  potentials.  The  values 
for  the  resistances  In  the  simple  model  are  determined  from  the 
values  for  which  * q^. 
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The  critical  parameters  In  the  determination  of  the  error 
due  to  axial  heat  flow  down  the  probe  are  the  probe  to  tissue 
resistance,  the  length  and  the  diameter  of  the  probe.  Blackwell 
(65,  66)  has  made  calculations  of  the  minimum  length  to  diameter 
ratio  for  thermal  probes  for  negligible  errors  due  to  axial  conduc- 
tion. He  calculates  steady  state  errors  for  the  axial  flow  problem 
as  a departure  from  the  solution  for  radial  flow  only.  He  concludes 
that  for  length  to  diameter  ratios  greater  than  25:1,  the  error  Is 
on  the  order  of  2 % for  solid  metal  probes  and  low  heat  transfer 
coefficient  values.  For  the  much  lower  conductivity  of  the  probes 
In  this  research,  the  error  should  be  much  smaller.  The  probe  In 
the  temperature  measurement  application  Is  In  1 millimeter  of  tissue 
and  the  diameter  is  less  than  20ym.  The  length  to  diameter  ratio 
Is  therefore  approximately  50:1,  and  the  steady-state  error  Is 
negligible. 

A final  property  of  the  probe  which  causes  errors  In  the 
measurement  of  transient  temperature  Increases  the  eye  or  any 
other  bio- system  Is  the  rise  time  of  the  probe  tissue  response.  j 

This  reponse  time  depends  upon  the  medium  In  which  the  probe  Is 
Imbedded,  the  thermal  properties  of  the  probe,  and  the  radius  of 
the  sensor  substrate.  The  resistance  between  the  probe  and  the 
tissue,  mentioned  In  the  previous  section,  also  affects  the  rise 
time.  From  measured  values  and  the  sizes  of  the  probes  used  In 
this  research,  the  response  time  for  the  probe  to  reach  90%  of  the 
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tissue  temperature  Is  less  than  one  millisecond.  This  value  is 
assumed  since  the  resistance  to  heat  flow  in  the  bio-system 
should  not  be  drastically  different  from  the  resistance  between 
the  probe  and  water  as  shown  in  Cain's  experimental  measurements 
(62).  For  10ym  radius  probes,  calculated  response  times  In  skin 
and  water  and  measured  response  times  are  all  less  than  one  milli- 
second for  Reed's  model  and  Cain's  experimental  measurements. 

Other  factors  In  the  system  which  limit  the  accuracy  of 
the  temperature  measurements  include  the  system  noise  level  and  the 
bandwidth  of  the  electronics  which  record  the  temperature.  The 


measured  input  noise  for  the  system  used  in  this  research  is 
4 yV RMS  referred  to  the  input  of  the  thermocouple  amplifier.  This 
gives  a resolution  of  ± .2°C  for  the  copper-nickle  thermocouple 
with  an  EMF  of  20yV/°C.  The  overall  bandwidth  of  the  measurement 
system  is  10  kHz,  which  is  equivalent  to  a rise  time  of  35yseconds 
for  a step  input. 

Other  characteristics  of  the  recording  system  also 
affect  the  accuracy  of  the  measurement  of  temperature  rise.  The 
analog  recorder  has  an  RMS  noise  level  of  .05V.  This  has  a 
varying  effect  on  the  temperature  measurement  since  the  gain  prior 
to  the  tape  system  is  Important  in  the  determination  of  the  minimum 
resolvable  temperature  change.  For  most  of  the  temperatures 
measured  In  the  results  section,  the  tape  noise  Is  less  than 
.5°C.  For  higher  values  of  pretape  gain,  the  error  is  much  less 
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than  .5°C.  The  analog  to  digital  conversion  adds  more  noise  to  the 
record  In  the  form  of  quantization  noise,  which  Is  ± 1/2  times  the 
least  significant  bit.  The  A/D  converter  is  a 12  bit  system.  The 
cumulative  effects  of  all  the  noise  In  the  recording  and  measurement 
system  amounts  to  less  than  ± 1°C.  This  represents  an  error  of 
± 5%  for  a 10  second  threshold  temperature  of  20°C. 


Threshold  Determination 


The  method  of  threshold  determination  in  this  research 
was  the  appearance  or  non-appearance  of  an  ophthalmoscoplcally 
visible  lesion  5 minutes  post  exposure  to  laser  Irradiation. 

This  criterion  was  chosen  to  allow  rapid  accumulation  of  data  and 
to  facilitate  threshold  determinations  for  several  exposure  dura- 
tions In  the  same  eye.  The  visible  lesion  criterion  is  relatively 
standard  throughout  past  research,  although  other  Investigators 
have  used  various  times  of  observation  post  exposure  as  the 
criterion  for  threshold  damage.  Even  though  the  appearance  of 
a visible  lesion  In  no  way  Implies  a loss  of  functional  vision. 

It  provides  a much  more  easily  determined  reference  point  for  the 
comparison  of  results  between  Investigators  than  histological  or 
hlstochemlcal  techniques.  For  the  reasons  of  easy  and  rapid 
threshold  determinations  and  the  availability  of  data  for  comparison, 
the  5 minute  visible  lesion  criterion  was  chosen. 
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Several  factors  affect  the  accuracy  of  threshold  deter- 
mination by  ophthalmoscopic  observation  techniques.  The  determina- 
tion of  the  visible  threshold  is  a subjective  technique  which 
depends  upon  the  condition  of  the  experimental  animal,  the  exper- 
tise of  the  observer  and  the  size  of  the  retinal  Image.  The  fundus 
of  the  experimental  animal  must  be  maintained  at  a reference 
temperature  (370C)  and  the  condition  of  the  eye  must  be  carefully 
monitored.  If  the  cornea  is  allowed  to  dry,  the  subsequent  clouding 
reduces  the  amount  of  light  reaching  the  retina  and  makes  visual 
observation  of  the  lesion  more  difficult.  In  this  research,  the 
temperature  of  the  fundus  is  kept  at  37°C  and  the  cornea  is  covered 
with  a contact  Tens  to  prevent  drying  of  the  corneal  epithelium. 

The  contact  lens  causes  a 1.6%  reduction  in  the  power  entering  the 
eye  since  its  index  of  refraction  is  higher  than  the  index  of 
refraction  of  the  water  layers  at  the  cornea. 

The  observer  makes  a judgment  of  lesion  appearance  based 
upon  changes  in  the  reflectance  of  the  fundus  following  laser 
exposure.  For  the  5 minute  visible  criterion,  a "threshold" 
lesion  is  one  which  appears  5 minutes  following  the  exposure. 

For  super- threshold  exposures,  the  lesion  may  appear  Immediately 
or  a short  time  following  the  Irradiation.  For  sub-threshold 
exposures,  lesion  appearance  may  or  may  not  occur  at  some  longer 
observation  time.  As  mentioned  In  the  background  section,  Davis 
and  Mautner  (14)  note  a 20%  difference  between  Immediately  visible 


threshold  lesions  and  lesion  appearing  one  hour  post  exposure. 
Gibbons  (15)  observes  no  difference  between  one  hour  and  24  hour 
lesion  criteria  for  times  of  exposure  less  than  5 seconds. 
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In  past  experimental  determinations,  the  corneal  power 
of  the  Irradlance  (J/cm2)  at  the  retina  have  been  used  as  a 
measurement  of  threshold  exposure.  The  range  of  power  or  energy 
for  a given  exposure  duration  varies  by  more  than  an  order  of 
magnitude  between  Investigators  and  within  a set  of  experiments 
by  a single  Investigator.  The  major  reason  for  this  discrepancy 
Is  probably  In  the  determination  of  retinal  Image  size.  The 
point  spread  function  of  different  animals,  even  of  the  same 
species,  exhibits  some  variation  since  not  all  animals  have  the 
same  total  refractive  power.  Residual  accoimodation  and  errors 
in  the  placement  of  external  lenses  can  cause  large  errors  in  the 
size  of  the  retinal  Image.  Errors  in  the  estimation  of  Image  size 
are  magnified  since  the  energy  per  unit  area  used  In  past  research 
Is  a function  of  the  square  of  the  assumed  radius  of  the  retinal 
Image.  In  the  measurement  of  thresholds  using  minimal  Images,  an 
error  of  ± .25  diopter  In  the  measurement  of  refractive  power  of 
the  eye  could  result  In  an  Image  of  34-45ym  rather  than  the 
assumed  10-20vim  Image  used  In  many  studies  (63).  For  this  reason, 
the  muscles  of  the  lens  must  be  paralyzed  and  an  accurate  measure- 
ment of  the  Image  size  must  be  made. 
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In  this  study,  paralysis  of  the  muscles  Is  accomplished 
with  curare  and  the  retinal  Image  size  Is  measured  directly  at  the 
retina  with  the  thermocouple.  The  effect  of  Image  size  and  distri- 
butlon  may  be  seen  from  the  tabulated  data  In  the  results  section. 
For  small  images,  the  power  required  Is  smaller  than  the  power 

required  for  large  Images.  For  Images  which  differ  significantly 
from  the  Gaussian  distribution,  the  power  required  to  produce  a 

lesion  changes  relative  to  a Gaussian  profile  with  the  same  half- 
power radius.  For  example,  a profile  with  large  tails  requires 
more  power  for  a threshold  exposure  than  a profile  with  a smaller 
amount  of  energy  in  the  tails  of  the  distribution. 

The  measurement  of  threshold  temperature  depends  upon  the 
linearity  of  temperature  rise  with  corneal  power  for  a given  Image 
distribution.  In  this  research  and  the  study  by  Cain  and  Welch 
(32),  the  linearity  of  temperature  versus  corneal  power  for  a 
given  image  distribution  has  been  demonstrated  for  temperature 
increases  up  to  50°C  above  baseline  temperature.  The  threshold 
temperature  in  most  of  these  experiments  is  extrapolated  from 
temperature  Increases  for  corneal  power  levels  50-75*  of  the 
threshold  value.  The  standard  deviation  of  the  threshold  tempera- 
tures for  each  exposure  duration  Is  within  15%  of  peak  temperature 
with  the  exception  of  .01  second  exposures  for  which  only  a few 
points  are  available.  The  corneal  power  standard  deviation  Is 
larger  than  the  standard  deviation  In  temperature  since  the  Image 
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size  and  distribution  appear  to  have  a much  greater  effect  on 
the  power  level  for  threshold  than  on  the  threshold  temperature. 

Lesion  Size  and  Lesion  Temperature  Measurements 

The  size  of  a lesion  which  Is  produced  by  exposure  to  an 
Intense  light  source  depends  upon  the  Image  size  and  the  duration 
of  the  exposure.  Also,  biological  amplification  may  play  a part 
In  the  lesion  size.  The  lesion  size  Is  larger  for  large  Images 
and  long  exposure  durations  than  would  be  expected  for  the  lesions 
produced  by  small  Images  and  short  duration  exposures.  Biological 
amplification  Is  due  to  changes  In  the  tissue  such  as  the  buildup 
of  toxic  metabolic  by-products  from  Increased  cellular  metabolism 
and  the  Immune  response  of  the  tissue  which  results  following 
cell  death.  The  lesion  may  therefore  be  larger  as  a result  of  this 
"amplification"  than  the  size  that  would  be  predicted  on  the  basis 
of  increased  temperature  alone.  In  fact.  In  response  to  Injury,  a 
minimum  volume  effect  due  to  cellular  necrosis  may  be  evident  for 
small  Images  or  short  exposure  durations  (63). 

Ophthalmoscopic  measurement  of  lesions  for  .01  second  to 
10  second  exposure  durations  are  listed  In  the  results  section  In 
Table  VI.  At  short  times  of  exposure,  the  radius  of  the  Image 
normalized  to  the  half-power  radius  Is  much  smaller  than  for  long 
times  of  exposure.  The  half-power  radius  Is  a better  parameter 
than  the  maximum  Image  radius  since  calculations  using  the  half-power 
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radius  as  the  reference  for  lesion  radius  result  In  much  lower 
standard  deviation  of  normalized  Image  radius.  Using  the 
ophthalmoscope  reticle  to  estimate  lesion  diameter  Is  limited  In 
accuracy  and  may  not  provide  a true  picture  of  the  lesion  behavior 
with  exposure  duration.  The  resolution  of  lesion  diameters  mea- 
sured with  the  ophthalmoscope  Is  at  best  ± 20ym.  The  resulting 
80-90%  accuracy  Is  due  to  the  13pm  separation  of  the  reticle 
divisions  and  the  lesion  boundaries  are  not  clearly  defined. 

The  lesion  appears  as  a gray  or  white  circular  area  which  repre- 
sents a change  In  the  reflectance  of  some  of  the  layers  In  the 
fundus.  The  lesion  probably  occurs  near  the  receptor  layers  (63) 
even  through  the  major  site  of  absorption  Is  In  the  pigment 
epithelium.  For  a lesion  to  be  visible,  coagulation  of  proteins 
probably  occurs  over  a large  volume  of  the  retinal  tissue  in  order 
to  produce  the  changes  In  the  diffuse  reflection  of  light  which 
make  the  observation  possible. 

The  fact  that  most  threshold  visible  lesions  require  some 
finite  time  to  develop  (at  least  for  exposure  durations  longer  than 
1 millisecond)  points  to  a biological  process  such  as  protein 
coagulation  and  denaturatlon  or  localized  edema  rather  than  a 
physical  response  such  as  steam  formation.  Around  the  melanin 
granules,  however,  the  temperature  may  be  locally  high  enough  to 
produce  steam.  This  Is  not  In  Itself  sufficient  to  produce  the 
effect  seen  In  the  formation  of  a visible  lesion,  since  most 
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lesions  are  not  Immediately  observable.  Localized  edema  and  fluid 
leakage  around  damaged  cells  certainly  causes  swelling  of  the 
tissue  and  a slight  change  in  the  refractive  index  of  the  tissue. 

All  these  effects  combine  to  produce  a visible  lesion  which  can 
be  observed  with  an  ophthalmoscope. 

Even  for  small  image  diameters,  a critical  volume  of 
tissue  must  be  involved  in  the  thermal  Insult  for  the  damage  to 
be  observable  with  an  ophthalmoscope.  The  observed  lesion  size 
for  minimal  Images  of  20pm  half-power  radius  Is  about  30-50pm. 

For  larger  images  the  size  of  the  damaged  volune  Increases  in 
radial  extent  (63).  The  minimum  lesion  size  for  small  Images  may 
be  a function  of  blolo^cal  amplification  or  a minimal  response 
to  tissue  Injury.  The  time  course  of  biological  amplification  due 
to  cellular  death  and  the  subsequent  buildup  of  toxic  byproducts 
follows  a time  course  similar  to  that  noted  for  lesion  formation 
(63).  The  damage  to  cells  In  the  area  of  a lesion  may  be  due  to 
a buildup  of  toxic  metabolic  byproducts  because  of  Interruption  of 
local  circulation  and  the  loss  of  nutrients  needed  to  sustain  the 
cell.  Thus,  cell  death  and  edema  probably  contribute  more  to  the 
appearance  of  a threshold  lesion  than  steam  formation  or  acoustic 
transients  In  the  range  of  exposure  durations  used  in  this  research. 
These  processes  are  the  result  of  elevated  temperature  and  depend 
upon  the  duration  and  Intensity  of  the  exposure  and  therefore  fit 
the  definition  of  a rate  process. 


Measurements  of  the  lesion  radius  are  used  In  conjunction 
with  the  horizontal  profiles  to  predict  the  temperature  rise  at 
the  lesion  radius.  The  lesion  radius  temperature  Is  extrapolated 
from  subthreshold  measurements  by  the  ratio  of  threshold  corneal 
power  to  the  corneal  power  for  the  particular  horizontal  scan. 

The  temperatures  at  a given  radius  for  a particular  Image  distribu- 
tion are  linear  with  corneal  power  as  are  the  center  temperatures 
discussed  earlier.  Because  of  the  inaccuracy  of  lesion  radius 
measurements,  the  variation  in  threshold  lesion  radius  temperatures 
is  larger  than  the  variation  In  lesion  center  temperatures.  The 
variation  In  lesion  radius  temperature  due  to  the  measurement  un- 
certainty is  much  smaller  for  10  second  temperature  measurements 
than  for  .01  second  temperatures  since  the  radial  gradients  are 
smaller  for  long  times  of  exposure  than  for  short  exposure  durations. 
The  uncertainty  in  lesion  radius  threshold  temperature  therefore 
becomes  less  as  the  duration  of  the  exposure  for  which  the  threshold 
Is  determined  Increases. 

Model  and  Experimental  Threshold  Comparisons 

In  the  comparison  of  model  and  experimental  threshold 
values,  a number  observations  may  be  made.  First,  real  differences 
exist  between  macular  and  paramacular  thresholds  determined  In 
this  study.  The  most  striking  difference  Is  the  difference  In 
temperature  Increase  required  for  visible  lesion  appearance  at  the 


two  Insertion  sites.  Another  difference  is  the  longer  rise  time 
for  paramacular  temperature  increases  compared  to  macular  Increases 
in  temperature.  The  model  calculated  threshold  temperatures  fit  the 
macular  experimental  values  more  closely  than  they  fit  the  para- 
macular experimental  threshold  temperature  increases. 

Since  the  damage  to  the  fundus  Is  a function  of  the 
temperature-time  history,  a longer  rise  time  for  the  temperature 
for  a given  image  size  or  Insertion  site  should  imply  a higher  end 
point  temperature  than  a temperature-time  history  with  a shorter 
rise  time.  Even  though  the  macular  rise  time  is  shorter  by  a 
factor  of  three  than  the  paramacular  rise  time  for  50ym  half-power 
images,  this  in  no  way  explains  the  approximately  20*  higher 
temperatures  required  for  paramacular  visible  lesions  if  the  simple 
rate  process  model  applies  to  the  Injury  of  retinal  tissues.  The 
rate  process  model  based  upon  the  Arrhenius  equation  is  a very 
sensitive  function  of  temperature.  For  the  20*  Increase  In  center 
temperature,  the  lesion  in  the  paramacular  region  should  be 
much  larger  in  radial  extent  and  axial  extent  than  the  lesion  In 
the  macula.  The  paramacular  threshold  lesion  therefore  would  Involve 
a much  larger  volume  of  tissue  than  the  macular  threshold  lesion. 
Unfortunately,  only  macular  lesion  measurements  have  been  made  for 
a number  of  experiments.  No  differences  can  be  demonstrated  In 
macular  and  paramacular  lesion  diameters  because  only  a few 
paramacular  lesions  have  been  measured  and  no  direct  comparisons  of 
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macular  and  paramacular  lesion  sizes  for  the  same  Image  size  have 
been  made. 

Several  observations  from  past  research,  however,  deserve 
further  note  In  regard  to  the  possible  reasons  for  macular  and 
paramacular  threshold  differences.  The  first  observations  relate 
to  structural  and  anatomical  differences  between  macular  and  para- 
macular regions.  The  macula  appears  darker  on  ophthalmoscopic 
observation  than  the  paramacula.  The  macula  also  contains  a yellow 
pigment  (Macula  Lutea)  which  absorbs  in  the  blue  and  green  spectra. 
Also,  the  neural  layers  are  much  thinner  In  the  area  of  the  macula 
and  the  fovea  than  in  the  paramacular  region  of  the  eye.  The 
dominant  cell  types  In  the  macula  and  paramacula  are  also  different. 
In  the  macula,  the  predominant  receptor  cell  is  the  cone  cell  while 
the  extramacular  receptor  cells  are  predominantly  rod  cells.  There 
is  more  circulation  in  the  neural  layers  of  the  paramacula  than  in 


the  macula.  The  neural  retina  In  the  paramacular  region  has  major 
vessels  and  extensive  capillary  circulation  while  the  macular  retina 
Is  largely  avascular.  Therefore,  several  major  anatomical  differ- 
ences exist  between  the  macula  and  paramacula.  Whether  these 
differences  contribute  to  the  observed  threshold  temperature 


rise  discrepancies  or  not  can  only  be  a matter  for  speculation. 

The  threshold  corneal  power  required  for  macular  and 
paramacular  lesions  produced  by  the  same  Image  sizes  Is  also  lower 
for  macular  thresholds  than  for  paramacular  thresholds.  For  50um 
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half-power  images  In  the  macula,  the  corneal  power  Is  approximately 
33  milliwatts  while  the  paramacular  threshold  is  65. 2 milliwatts 
for  a .1  second  exposure  duration.  For  a single  threshold  in  the 
same  eye,  with  105pm  half-power  radius,  the  macular  lesion  requires 
a corneal  power  of  72.1  milliwatts  while  the  paramacular  lesion 
requires  101  milliwatts  for  threshold  corneal  power  at  .1  second 
exposure  duration.  At  one  second,  macular  lesions  require  20-25% 
less  corneal  power  than  paramacular  lesions.  Ten  second  macular 
threshold  lesions  are  produced  with  20-25%  less  corneal  power  than 
paramacular  threshold  lesions.  These  observations  parallel  the 
results  of  other  investigators  (6,  12,  67)  for  other  wavelengths. 
Lappin  (12)  found  a reduction  in  corneal  power  for  macular  lesions, 
for  He-Ne  wavelengths,  at  1 to  10  seconds  of  20-25%  from  the  power 
required  for  threshold  lesions  at  extramacular  exposure  sits.  At 
.01  second,  the  paramacular  lesion  threshold  corneal  power  is  more 
than  twice  the  macular  threshold  power.  For  ruby  exposures, 
Vassilladls,  et  al.  (6)  find  that  the  macula  Is  2.2  times  more 
sensitive  than  the  paramacula  for  1.7  millisecond  exposures.  They 
state  that  the  reason  for  the  difference  may  be  that  the  macula 
Is  more  heavily  pigmented  than  the  extramacula.  This  Implies  that 
the  absorption  of  light  in  the  macular  P.E.  and  choroid  is  greater 
than  the  absorption  In  the  paramacula.  Campbell  (67)  finds  the 
macula  2.5  times  more  sensitive  to  ruby  laser  Irradiation  than  the 
paramacular  for  exposure  durations  shorter  than  1 millisecond. 
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The  reasons  for  the  difference  In  sensitivity  of  the 
macular  and  extramacular  regions  of  the  eye  Is  not  clear,  but 
generally  the  arguments  for  the  difference  tend  to  support  several 
hypotheses:  (1)  That  the  observation  of  a lesion  In  the  macula  Is 

easier  than  the  observation  In  the  paramacula  due  to  the  darker 
background,  (2)  That  the  absorption  profile  In  the  macula  and  para- 
macular regions  are  different,  (3)  The  cones  and  rods  may  exhibit 
differential  sensitivity  to  damage  by  Intense  light  which  Implies  a 
different  damage  mechanism  for  macular  and  paramacular  regions,  and 
(4)  anatomical  differences  between  the  macula  and  paramacula  con- 
tribute to  the  different  threshold  temperatures  In  the  two  areas. 

The  differences  In  macular  and  paramacular  light  absorption 

and  thickness  probably  contribute  to  the  darker  appearance  of  the 

macula.  Brindley  (68)  finds  a difference  in  optical  density  of  the 

o 

neural  retina  In  the  two  regions  of  a factor  of  two  at  5000A  wave- 
lengths. The  macular  retina,  even  with  Its  reduced  thickness  has  a 
density  twice  that  of  the  paramacular  neural  retina.  This  factor 
could  account  for  the  darker  appearance  of  the  macula  ophthalmo- 
scoplcally.  The  change  In  reflectance  which  represents  a lesion 
could  therefore  be  a result  of  the  requirement  that  a larger 
volume  of  tissue  must  be  denatured  or  be  edematous  In  the  thicker, 
less  absorbent  extramacular  retina  than  In  the  thinner,  more 
absorbent  macular  neural  retina.  A smaller  volume  of  tissue 
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change  would  produce  the  same  reflection  In  the  macula  that 
a larger  volume  of  tissue  change  would  produce  In  the  paramacula 
for  wavelengths  near  500oR.  The  thicker  paramacular  retina  causes 
more  scattering  of  light  since  for  a small  volume  more  retinal 
thickness  is  present  between  the  observer  and  the  lesion  boundary. 
Assuming  that  less  light  is  transmitted  through  the  layer  from  the 
lesion  in  the  paramacula  because  of  light  scatter  than  in  the 
thinner  macular  layers,  a small  lesion  or  tissue  volume  change 
would  be  more  observable  in  the  macula.  For  longer  wavelengths, 
the  two  areas  exhibit  approximately  the  same  optical  density. 

The  second  hypothesis  of  higher  absorption  in  the  P.E. 
and  choroid  of  the  macula  than  in  the  extramacular  region,  may  be 
supported  by  the  observation  that  higher  corneal  power  is  required 
to  produce  a threshold  lesion  in  the  paramacular  areas  of  the  eye. 
This  assumption  does  not  Indicate  that  the  temperature  should  be 
different  for  the  two  sites,  but  the  absorption  parameters  must  be 
different.  For  short  exposure  durations,  the  heat  conduction 
equation  becomes: 

qHP  . 3T 

pc  3 1 

and  the  temperature  rise  Is  proportional  to  the  absorption  parameter 
and  the  Irradlance  HQ.  The  longer  rise  time  of  the  paramacular 
region  Indicates  that  it  has  a smaller  absorption  coefficient  than 
the  macular  region  of  the  eye. 
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The  difference  In  corneal  power  required  for  macular  and 
paramacular  lesions  appears  to  be  present  at  argon  wavelengths 
used  In  this  research  as  well  as  the  helium  neon  (12)  and  ruby  (6) 
wavelengths.  This  Indicates  that  some  difference  exists  In  either 
the  volume  of  damage  required  to  produce  an  observable  effect  In 
the  two  major  areas  of  the  eye  or  of  some  difference  In  the  damage 
mechanism  by  which  the  lesion  becomes  visible.  Tso  et  al.  (69, 

70,  71)  In  studies  of  the  repair  and  development  of  retinal  lesions 
find  that  cone  cells  are  much  more  susceptible  to  damage  than  rod 
cells.  Although  their  studies  are  performed  with  an  argon  laser 
for  the  evaluation  of  differential  rod  and  cone  damage  and  they 
are  done  in  the  periphery  of  the  monkey  eye,  they  feel  that  the 
susceptibility  of  the  cone  cells  Is  not  wavelength  dependent. 

They  conclude  that  the  differences  In  susceptibility  of  the 
two  cell  types  to  laser  radiation  explain  the  Increased  corneal 
power  required  to  produce  lesions  In  macular  and  paramacular 


areas  (70). 

Since  experiments  of  Tso  et  al.  are  conducted  In  areas 
of  the  retina  which  contain  both  rod  and  cone  receptors,  the  ob- 
served differences  In  damage  to  the  two  cell  types  are  probably 
due  to  general  morphological  differences  noted  In  their  histologic 
studies.  First,  cone  and  rod  cells  exhibit  different  Interfaces 
with  the  pigment  epithelium.  The  cone  cells  are  surrounded  by 
more  densely  plgnented  villi  from  the  P.E.  than  are  the  rods. 


The  villi  from  the  P.E.  also  appear  to  extend  farther  around  the 
cone  cell  outer  segments  than  around  the  rod  outer  segments.  Nuclei 
of  cone  cells  are  observed  often  in  the  inner  segments  of  the 
receptor  layer  in  the  normal  retina  as  an  anatomic  variant. 

The  nuclei  of  cone  cells  also  appeared  to  be  more  sensitive  to 
thermal  effects  from  photic  injury  than  the  nuclei  of  rod  cells 
since  the  nuclei  of  the  cones  are  physically  nearer  to  the  pigment 
epithelium.  The  outer  segments  of  cones  In  detached  retinas  of 
rhesus  monkeys  are  observed  to  recover  more  slowly  than  the  outer 
segments  of  rods. 

The  lesions  in  their  studies  extend  over  the  inner 
and  outer  segments  of  the  receptor  layer  of  the  neural  retina  with 
the  approximate  center  at  the  external  limiting  membrane.  All 
nuclei  in  the  outer  nuclear  layer  are  affected  initially,  and  the 
swelling  of  the  tissue  due  to  edema  causes  normal  nuclei  to  be 
displaced  from  the  lesion  center.  In  the  periphery  of  the  lesion 
only  rod  nuclei  may  be  seen  at  long  times  post  exposure.  The  rod 
cells  apparently  regenerate  after  some  period  of  time  with  the 
cone  cells  being  selectively  destroyed.  This  effect  Is  also  seen 
In  Xenon  arc  lamp  coagulation  which  leads  to  the  conclusion  that 
the  sensitivity  of  rod  cells  may  not  be  wavelength  dependent. 

In  the  macular  region  of  the  eye,  the  cone  cells  are 
densely  packed  and  no  neural  layer  circulation  exists.  In  the 
fovea,  an  extremely  high  cell  density  of  cones  is  present.  As 
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the  distance  from  the  macula  Increases,  the  density  of  cones  Is 
very  low  and  rods  are  predominant  cell  types.  In  the  periphery 
of  the  retina,  rod  cells  of  lower  density  are  predominant.  In 
peripheral  and  extramacular  neural  layers,  the  Inner  layers  have 
an  extensive  capillary  circulation. 

The  absorptive  properties  of  blood  and  blood  vessels  are 
higher  for  argon  wavelengths  than  for  other  visible  wavelengths. 

Thus,  In  the  paramacula,  some  of  the  energy  at  argon  wavelengths 
is  lost  In  the  retinal  layers.  However,  as  even  larger  percent 
of  light  may  be  lost  in  the  Macula  Lutea  of  the  macular  region. 

These  source  terms  are  not  accounted  for  in  the  model  since  no 
reliable  measurements  of  the  magnitude  of  this  absorption  are 
available. 

The  difference  in  threshold  power  for  macular  and  para- 
macular regions  may  be  explained  by  damage  volume  differences 
required  for  observations  or  by  the  differential  absorption  and 
sensitivity  of  cone  and  rod  areas  of  the  retina.  It  does  not 
appear  that  a true  explanation  of  the  differences  In  temperature 
rise  can  be  resolved  until  histological  or  hlstochemlcal  evidence 
of  damage  volume  differences  Is  available.  The  apparent  differential 
absorption  phenomena  also  need  to  be  Investigated  for  both  macular 
and  paramacular  areas.  Past  studies  have  not  Indicated,  In  suf- 
ficient detail,  adequate  reasons  for  the  temperature  difference 
noted  In  this  research.  Without  histological  evidence  and  accurate 


measurements  of  absorption  parameters,  the  model  cannot  be  expected 
to  describe  the  effects  that  have  been  noted  in  this  and  other 
research.  It  appears,  however,  that  because  of  real  anatomical 
differences  between  the  two  major  areas  of  the  eye  that  a simple, 
universal  model  for  temperature  Increase  and  thermal  damage  cannot 
adequately  describe  the  Injury  In  the  two  areas. 

Model  and  Experimental  Scans 

From  the  horizontal  and  axial  scan  data  compared  with 
model  values,  several  differences  between  model  and  experimental 
values  are  apparent.  First,  horizontal  scans  in  diacular  and 
paramacular  areas  reflect  the  problems  noted  In  the  comparison  of 
temperature-time  histories.  In  the  Results  section,  the  measured 
rise  time  for  the  paramacular  temperature  rise  is  seen  to  be 
longer  than  the  rise  time  for  temperature  Increases  in  the  macular 
region. 

Generally,  the  paramacular  measured  center  temperature 
are  higher  than  the  model  predictions  for  1 second  and  10  second 
exposure  durations.  For  exposure  durations  of  .1  second  and  below, 
the  measured  temperature  Increases  are  lower  than  the  calculated 
values.  However,  at  large  radial  displacements  from  the  lesion 
center,  the  measured  paramacular  temperatures  are  higher  than  the 
model  temperatures.  This  effect  Is  normally  at  displacements  much 
larger  than  the  lesion  radius.  It  may  be  caused  by  the  presence 
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of  blood  flow  In  the  paramacular  region  which  causes  the  effective 
conductivity  to  be  larger  since  blood  which  absorbs  light  more 
heavily  In  the  argon  spectrum  than  In  other  parts  of  the  spectrum 
enhances  the  transport  of  heat  by  a convective  mechanism.  In  the 
macular  region  the  lack  of  blood  flow  Is  probably  the  reason  that 
this  effect  is  not  noted  for  horizontal  scans  in  the  macula. 

Axial  scans  in  the  two  areas  of  the  eye  are  much  smoother 
than  the  model  predicted  scan  for  parameter  set  1.  Some  of  this 
effect  may  be  caused  by  probe  averaging  as  mentioned  earlier. 

The  axial  scans  are  the  most  difficult  of  the  scans  to  accomplish. 
Some  of  the  problems  which  may  cause  discrepancies  between  model 
and  experimental  axial  scans  are:  (1)  the  condition  of  the  probe 

track  through  the  tissue;  (2)  the  probe  may  hang  on  tissue  during 
the  scan;  (3)  the  assumed  homogeneity  versus  the  actual  Inhomo- 
genous  structure  of  the  various  layers  in  the  fundus,  and  (4)  the 
lack  of  accurate  knowledge  of  the  absorption  parameters  In  the 
various  layers  and  areas  of  the  fundus.  The  first  point  is  very 
Important  to  any  attempt  to  measure  the  absorption  parameters  In 
the  tissue  by  means  of  an  axial  scan  and  to  the  measurement  of 
temperature  In  the  various  layers.  The  track  of  the  probe  through 
the  tissue  Is  at  an  angle  to  the  axis  of  the  beam  and  the  finite 
size  of  the  probe  causes  some  tissue  deformation.  If  a hole  In 
the  tissue  exists,  a larger  than  normal  direct  absorption  component 
will  be  present  on  the  axial  temperature  Increase  for  at  least  the 
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first  several  tens  of  micrometers  behind  the  P.E.  A subretlnal 
hemmorhage  In  the  choroid  will  cause  a lower  blood  perfusion  rate 
than  for  normal  tissue  and  may  cause  a local  change  In  the 
absorption  of  light.  A hole  In  the  neural  layers  or  In  the  P.E. 
and  choroid  Increases  the  contact  resistance  of  the  probe-tissue 
Interface.  Host  of  these  effects  probably  cause  perturbation  of 
temperature  fields  In  the  tissue,  especially  for  short  exposure 
durations.  The  volume  effect  of  tissue  heating  tends  to  reduce 
these  problems  for  longer  exposure  durations  If  the  abnormality 
In  the  tissue  Is  confined  to  a small  volune  of  the  heat  source. 

In  most  scans,  the  probe  Is  withdrawn  from  the  neural 
layers  in  small  steps  rather  than  driven  Into  the  tissue  from  the 
sclera.  The  reasons  for  this  approach  are  that  the  Image  must 
first  be  located  and  the  simplest  method  Is  by  visual  positioning 
of  the  probe  in  the  laser  Image  through  the  fundus  camera  and 
that  the  probe  hangs  on  the  tissue  more  frequently  in  the  forward 
direction.  If  the  probe  hangs  In  the  reverse  direction,  the  scan 
Is  not  good,  but  If  It  hangs  In  the  forward  direction  it  usually 
results  In  a broken  probe.  If  the  probe  hangs  on  the  tissue  during 
withdrawal,  the  temperature  exhibits  an  abrupt  decrease  at  some 
point  during  the  scan. 

The  model  assumes  that  the  tissue  Is  homogeneous,  with 
constant  absorption  parameters  In  the  P.E.  and  choroid  layers. 

This  Is  clearly  not  an  accurate  representation  of  the  tissue  system. 


The  tissue  Is  not  arranged  In  homogeneous  layers.  In  fact, 
uncertainty  In  the  definition  of  the  depth  of  the  layers  as  well 
as  local  variations  In  the  blood  flow,  pigmentation,  and  layer 
thickness  make  the  descriptions  of  the  absorption  profiles  and 
therefore  the  variation  In  temperature  versus  position  for  short 
exposures  very  difficult.  The  P.E.  is  not  a homogeneous  10pm  thick 
layer  as  the  model  assumes.  This  layer  Is  non-uniform  In  thickness 
and  has  projections  (villi)  which  reach  into  the  neural  layers 
(outer  segments)  which  vary  In  pigmentation  and  length  around  rod 
and  cone  receptors  as  discussed  before  (70).  The  choroid  layer 
contains  both  blood  vessels  and  pigment  granules  which  vary  as  a 
function  of  position  In  the  eye.  The  choriocapillarls  which 
underlies  the  P.E.  extends  only  several  tens  of  microns  Into  the 
lGOym  assumed  thickness  of  the  choroid.  The  band  of  pigmentation 
Is  Irregular  in  extent  through  the  choroid  layer.  These  observa- 
tions Indicate  that  more  layers  should  be  Included  in  the  model 
especially  for  short  duration  exposures  for  which  energy  density 
rather  than  volume  effects  are  Important.  The  assumption  of 
homogeneity  may  not  be  too  bad  for  a point  near  the  center  of  the 
Irradiated  volume  of  tissue  for  long  exposure  times.  The  exposure 
durations  used  In  this  research  are  relatively  long  compared  to 
the  times  required  for  significant  conduction  effects  to  "average" 
the  affects  of  varying  absorption  within  layers.  For  times  longer 
than  1 millisecond  the  tissue  can  probably  be  represented  by  average 


absorption  and  thermal  parameters  within  some  small  range  of 
error. 

Model  Parameters 

The  value  of  the  parameters  In  the  thermal  model  are 
very  Important  to  the  accurate  prediction  of  temperature  Increase 
in  the  fundus.  As  noted  before,  the  absorption  profile  In  the 
tissue  Is  probably  not  the  same  In  the  tissue  as  that  assumed  for 
homogeneous  absorbing  tissue  layers.  For  long  exposures,  however, 
the  averaging  effect  of  heat  conduction  tends  to  reduce  the  effect 
of  small  Inhomogeneities.  At  short  times,  the  Important  parameters 
are  absorption  coefficients  and  volumetric  specific  heats  for  the 
absorbing  bodies  in  the  tissue.  Calculations  by  Hayes  and 
Wolbarsht  (72)  indicate  that  the  melanin  granules  may  achieve 
temperature  increases  of  several  hundred  degrees  centigrade  for 
theshold  exposures  below  a few  microseconds.  If  this  Is  true,  local 
variations  In  tissue  absorption  around  the  melanin  granules  surely 
occur  due  to  denaturatlon  and  coagulation  of  protein  and  steam 
bubble  production  In  the  tissue  surrounding  the  granules.  The 
.01  second  duration  exposures  of  tf.<s  research  may  exceed  phase 
change  temperatures  near  the  center  of  the  lesion.  Change  of  phase 
may  also  cause  changes  In  the  tissue  absorption  properties.  Thus, 
the  assumption  of  constant  absorption  parameters  for  the  various 
layers  Is  probably  not  valid  for  short  duration  threshold 
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temperatures  with  the  model. 

Local  blood  flow  effects  In  the  macula  and  paramacula 
are  not  accounted  for  In  the  model.  The  effect  of  blood  flow  on 
the  temperature  Increase  Is  not  even  Included  In  the  slnpllfled 
thermal  model  used  to  predict  retinal  temperature  Increase.  The 
parameters  Involved  In  the  analysis  of  blood  flow  should  account 
for  the  Influence  of  not  only  the  cooling  effect  on  the  tissue 
but  also  the  spectral  absorption  of  the  vessels  and  hemoglobin 
within  the  layers  In  which  significant  blood  flow  Is  present. 

No  accurate  numbers  exist  for  either  blood  perfusion  rates  or  the 
absolute  value  of  hemoglobin  spectral  absorption.  Blood  flow  Is 
variable  within  the  various  regions  and  layers  of  the  tissue. 

Beble  et  al.  (73)  estimate  that  22%  of  the  argon  laser  light 
reaching  the  paramacular  neural  retina  Is  absorbed  In  the 
capillary  bed  in  the  neural  layers. 

Measurement  of  direct  absorption  on  the  thermocouple  as 
a function  of  axial  position  in  the  tissue  does,  however,  appear 
to  fit  the  exponential  absorption  profile  in  the  results  section 
relatively  well.  The  probe  smoothing  effect  and  probe  track  condi- 
tion could,  however,  have  a great  effect  on  these  measurements. 

The  calculated  absorption  coefficients  for  two  scans  are  648  cm-1 
and  119.4  cm-1  for  the  macular  P.E.  and  choroid;  for  one  run  In  the 
paramacula  P.E.  and  choroid  coefficients  appear  to  be  390.2  cm'1 
and  102.6  cm"1.  The  small  number  of  scans  are  due  to  the  spectral 
techniques  which  must  be  employed  to  obtain  these  scans.  The  probe 
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track  must  be  In  good  condition*  the  probe  cannot  hang  on  the 
tissue,  and  a high  power  and  short  pulse  must  be  used  to  give 
sufficient  direct  absorption  for  the  measurement. 

The  measured  parameters  dp,  however,  closely  correspond 
to  parameter  set  2 for  the  macular  experimental  values  and  to 
parameter  set  3 for  the  paramacular  experimental  parameters.  If 
these  parameters  are  accurate,  they  provide  a partial  explanation 
for  the  differences  In  average  rise  time  between  the  macular  and 
paramacular  temperature-time  histories.  The  higher  absorption 
parameter  Implies  a greater  initial  slope  of  temperature  versus 
time  for  the  macula  and,  therefore  a shorter  rise  time.  The  model 
also  reflects  this  for  parameter  sets  2 and  3.  Parameter  set  2 has 
a greater  Initial  slope  of  temperature  versus  time  and  a faster  rise 
time  due  to  the  higher  value  for  the  P.E.  absorption  coefficient. 

In  Table  XVI,  the  average  values  of  threshold  temperature 
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In  Table  XVII,  the  average  lesion  threshold  center  temperatures 
for  paramacular  exposures  are  compared  to  model  predicted  values 
for  10  second  to  .1  second  exposure  durations.  It  appears  that 
model  set  3 provides  a slightly  better  "fit"  to  the  macular 
experimental  data  and  a much  better  fit  for  the  paramacular  data. 

TABLE  XVII:  Comparison  of  Paramacular  Experimental 

and  Model  Threshold  Temperatures 

Exposure  Experimental  Model  Set  1 Model  Set  2 Model  Set  3 
Duration  Temperature  Temp.  Rise  Temp.  Rise  Temp.  Rise 

(seconds)  Rise  (°C) (°C) (°C) (°C) 

10  24.8  ±1.7  21.9  ± .6  19.9  ± 5.4  24.9  ± 5.3 

1 28.7  ± 3.0  27.3  ± 4.5  25.1  ± 5.8  29.4  ± 4.9 

.1  40.0  ± 5.1  47.7  ± 9.0  43.5  ± 9.6  45.5  ± 8.2 


Since  paramacular  exposures  represent  only  about  one-third  of  the 
exposures  In  this  research  direct  comparisons  are  difficult.  The 
average  temperature  Increases  do  not  reflect  changes  In  Image  sizes 
between  macular  and  paramacular  Insertions  as  the  corneal  power 
average  values  do.  The  relatively  good  agreement  of  the  model  for 
both  macular  and  paramacular  exposures  of  duration  longer  than  .1 
second,  for  all  parameter  sets.  Indicates  that  the  volume  of  the 
source  term  rather  than  the  absolute  value  of  the  absorption  para- 
meters, Is  most  Important  In  the  determination  of  tissue  temperature 


Increase. 


The  absorption  values  for  set  1 are  probably  too  high  for 
the  preretlnal  ocular  media  and  the  P.E.  Sets  2 and  3 have  a more 
realistic  value  for  TOM  based  upon  earlier  discussions.  The  values 
of  absorption  in  the  P.E.  and  choroid  for  sets  2 and  3 agree  fairly 
closely  with  the  values  obtained  for  a small  number  of  axial  scans 
in  the  macula  and  paramacula  respectively.  The  model  exhibits  good 
agreement  with  both  macular  and  paramacular  temperature  increases 
over  the  range  of  exposure  durations  from  .1  second  to  10  seconds. 
The  values  obtained  are  within  the  range  of  experimental  errors 
discussed  in  earlier  sections. 

At  .01  seconds,  several  problems  were  noted  during 
the  course  of  the  experiments  which  affect  the  threshold  temperature 
measurements  and  model  values.  First,  for  the  higher  corneal  power 
required  to  produce  a lesion  at  .01  seconds,  the  neutral  density 
filters  heated  causing  a change  in  power  with  time.  Second,  the 
laser  exhibited  larger  divergence  angles  for  high  power  than 
for  lower  power  settings.  These  problems  were  not  present  for  the 
lower  power  settings  required  for  durations  of  exposure  longer  than 
40  milliseconds.  The  .01  second  exposures  therefore  required  a 
different  optical  system  and  a much  more  time  consuming  experiment 
than  the  longer  duration  exposures.  The  large  standard  deviations 
seen  In  the  .01  second  exposures  were  due  to  the  errors  induced  by 
he  problems  mentioned  above. 
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The  model  of  Henri ques  (51)  does  not  accurately  describe 
the  damage  to  retinal  tissues.  The  coefficients  derived  for 
macular  and  paramacular  areas  are  different  from  one  another  and 
are  not  the  same  as  Henrlques  coefficients.  The  reason  for  these 
observed  differences  In  damage  temperatures  have  been  discussed  In 
a previous  section.  The  major  difference  Is  probably  due  to  the 
fact  that  the  damage  observed  In  this  study  Is  not  the  result  of 
a simple,  first  order  reaction.  The  damage  observed  Is  probably 
due  to  denaturatlon  and  coagulation  of  retinal  proteins  and  tissue 
ea«ir.a  which  results  from  damage  to  retinal  cells  and  deactivation 
of  enzyme  systems  which  are  important  In  the  synthesis  of  retinal 
proteins.  Of  these,  tissue  edema  is  probably  the  major  effect 
which  causes  the  observed  change  In  fundus  reflectance. 

The  rate  process  model  depends  upon  a knowledge  or  the 
ability  to  predict  temperature  at  all  positions  In  the  tissue. 

The  process  by  which  the  damage  occurs  must  be  known  or  modeled 
empirically  from  observations  of  damage  by  some  arbitrary  criterion. 
At  best,  because  of  experimental  error  and  variations  In  criteria 
for  damage  only  an  average  effect  may  be  predicted.  Especially, 
the  apparent  volume  effect  of  damage  rather  than  the  damage  at  a 
point  must  be  considered.  The  possibility  that  the  tissue  exhibits 
a minimum  response  to  tissue  damage  which  Is  a function  of  biological 


amplification  that  may  not  correspond  to  the  same  criteria  as  for 
protein  or  enzyme  changes  must  be  considered.  Errors  In  the  estima- 
tion of  damage  may  be  significant  If  a minimum  volume  must  be 
damaged  In  order  to  observe  a lesion.  If  the  system  response  (e.g. 
edema)  Is  larger  than  the  actual  volume  of  damage,  a rate  process 
model  based  on  temperature  Increase  alone  may  not  predict  the 
observed  lesion  radius  accurately. 

The  retinal  damage  produced  by  threshold  exposures  from 
.01  second  to  10  seconds  In  this  study  can  be  predicted  by  a rate 
process  model.  (1)  The  temperature  required  to  produce  the  effect 
which  Is  the  criterion  for  damage  Increases  with  decreasing  exposure 
duration.  (2)  The  radius  of  the  visible  lesion  decreases  with 
decreasing  exposure  time  due  to  the  reduced  effect  of  conduction 
at  shorter  exposure  durations.  Even  though  the  measurement  of 
lesion  radius  and,  therefore,  the  actual  temperature  at  which 
damage  occurs  cannot  be  measured  accurately  by  ophthalmoscope, 
the  trends  which  Imply  that  the  damage  observed  Is  a thermal  rate 
process  are  clearly  demonstrated  by  this  research. 
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CHAPTER  VI 
CONCLUSIONS 

Temperature  Increases  In  the  primate  fundus  (Hacaca 
mulatta)  for  threshold  lesion  appearance  5 minutes  post  exposure 
to  argon  laser  Irradiation  have  been  measured  for  exposure  durations 
from  .01  second  to  10  seconds.  Comparisons  of  measured  and  cal- 
culated temperature  Increases  have  been  presented  for  three  sets 
of  model  parameters  In  the  computer  solution  of  the  heat  conduction 
equation.  The  rate  process  model  predictions  of  damage  temperature 
have  been  compared  with  the  measured  lesion  center  and  lesion  radius 
temperature  Increases.  Conclusions  resulting  from  this  research 
are  as  follows: 

1.  The  measured  Image  center  temperature  Increases  at 
the  axial  position  of  highest  temperature  rise  Is 
higher  for  paramacular  exposures  than  for  macular 
exposures.  The  mean  and  standard  deviation  of 
temperature  Increase  for  paramacular  exposure 
durations  of  .01  second,  .1  second,  1 second,  and 
10  seconds  are  S5.9  ± 14.9°C,  39.9  i 5.1°C, 

28.7  ± 3.0°C,  and  24.8  ± 1.7°C,  respectively.  For 
macular  exposures  of  .01  second,  .1  second,  1 
second,  and  10  seconds  duration,  the  average 
values  and  the  standard  deviations  of  the  measured 
temperature  Increases  are  66.1  ± 44.3°C,  29.4  ± 4.1°C, 
180 
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24.1  ± 3.4°C  and  19.9  ± 1.5°C  respectively. 

2.  The  damage  temperatures  at  the  ophthalmoscoplcally 
measured  lesion  radius  are  estimated  to  be  .58,  .75, 
.82  and  .88  times  the  lesion  center  temperatures  for 
macular  and  paramacular  exposure  durations  of  .01 
second,  .1  second,  1 second  and  10  seconds  respective- 
ly. The  error  In  this  measurement  Is  lower  for  long 
duration  exposures  than  for  short  duration  exposures. 

3.  The  temperature  model  predicts  lesion  center  tempera- 
tures within  20*  of  the  measured  temperature  Increases 
for  exposure  durations  from  .1  second  to  10  seconds 
for  paramacular  exposures  and  within  10*  for  macular 
exposures.  The  model  with  parameter  set  3 is  within 
14*  for  paramacular  exposures  ar.d  6*  for  macular 
exposures  for  exposure  durations  of  .1  second  to  10 
seconds. 

4.  The  damage  to  retinal  tissues  resulting  from  laser 
Irradiation  for  exposure  durations  used  In  this 
research  Is  the  result  of  Increased  temperature  In 
the  retinal  tissues.  Exposure  durations  shorter  than 
10  seconds  required  a higher  temperature  Increase  to 
produce  damage.  Therefore,  the  damage  to  the  tissue 
conforms  to  the  definition  of  a rate  process  since 
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the  damage  Is  dependent  upon  both  the  temperature 
Increase  and  the  exposure  duration. 

5.  A simple  rate  process  for  protein  denaturatlon  by 
Henrlques  does  not  predict  the  estimated  damage 
temperature  at  the  lesion  radius  found  In  this 
research.  In  the  macula,  the  predicted  damage 
temperatures  are  higher,  while  In  the  paramacular, 
the  predicted  temperatures  are  lower  than  the 
experimental  values.  This  Implies  a volume  difference 
for  damage  which  results  In  an  observable  lesion  with 
the  damaged  volume  being  larger  In  the  paramacular 
region  of  the  eye. 

6.  Both  measured  and  modeled  temperatures  display  a 
lower  standard  deviation  than  the  corneal  power  stand- 
ard deviation.  This  Implies  that  the  damage  tempera- 
ture Is  less  sensitive  to  changes  In  Image  radius 

and  the  shape  of  the  Image  distribution  than  corneal 
power  over  the  limited  range  of  Image  sizes  and 
distribution  In  this  researchr'Therefore , an  accurate 
knowledge  of  Image  size  and  distribution  of  power  In 
the  Image  Is  necessary  for  comparison  of  corneal 
powers  required  to  produce  damage. 
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7.  The  difference  In  rise  time  between  the  model  and 
the  experiment  shows  a need  for  the  accurate  deter- 
mination of  the  absorption  parameters  for  the  P.E. 
and  choroid.  For  the  model  to  accurately  predict 
temperature- time  histories  and  damage  profile,  the 
time  constant  as  well  as  "end  point"  temperature 
must  be  accurately  matched  since  the  rate  process 
model  employs  an  exponential  function. 

8.  From  this  research.  It  appears  that  a universal 
model  (l.e.,  constant  absorption  coefficients) 
cannot  be  used  tc  accurately  predict  retinal  tempera- 
ture In  the  macular  and  paramacular  regions  of  the 
eye.  The  differences  In  both  the  rise  time  of 
temperature  and  the  final  value  which  produce  damage 
in  the  two  areas  imply  that  different  coefficients 
and  damage  models  are  required  for  the  prediction  of 
fundus  damage. 

9.  Further  experiments  are  necessary  to  provide  accurate 
predictions  of  temperature  rise  at  the  lesion  radius 
and  the  extent  damage  In  the  retina  which  results 
from  exposure  to  Intense  light  sources.  Some  of  the 
basic  questions  which  require  further  Investigation 
are: 
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(a)  The  reasons  for  differences  In  macular  and 
paramacular  threshold  temperature  Increases  must 
be  determined.  Histological  and  hlstochemlcal 
studies  are  needed  to  determine  macular  and 
paramacular  differences  at  the  ophthalmoscopic 
thresholds  for  a number  of  exposure  durations. 

In  vivo  measurements  of  the  absorption  of  light 
In  the  neural  layers,  pigment  epithelium,  and 
choroid  are  required  as  a function  of  wave- 
length and  position  In  the  fundus.  Accurate 
determinations  of  the  preretlnal  ocular  media 
transmission  and  the  absorption  coefficients  of 
the  various  layers  of  the  fundus  are  necessary 
for  the  accurate  prediction  of  temperature 
Increase  In  these  two  areas  of  the  eye. 

(b)  Determination  of  the  effect  of  blood  flow  on  the 
retinal  temperature  Increase.  Measurements 

of  threshold  temperatures  with  and  without  blood 
flow  can  be  used  to  determine  the  Influence  of 
vascular  perfusion  on  temperature  rise  time  and 
the  final  temperature  required  to  produce  damage. 
Accurate  estimates  of  blood  flow  rates  and  the 
effect  of  blood  flow  on  model  calculated  tempera- 
ture-time histories  are  also  needed. 
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The  optical  transfer  function  for  the  prediction 
of  damage  from  minimal  Image  sizes  must  be  ob- 
tained. The  model  requires  accurate  determina- 
tion of  the  point  spread  function  of  the  eye  and 
measurement  of  temperatures  for  minimal  Image 
lesions.  Histological  studies  should  also  be 
used  to  determine  the  effect  of  Image  size  on 
damage  volume  for  an  observable  lesion. 
Determination  of  changes  In  functional  vision 
as  evidenced  by  reduced  local  electrical  activity 
for  threshold  and  sub-threshold  lesions.  Local 
potential  measurements  In  the  lesion  from  a 
combination  thermocouple  and  potential  measurement 
electrode  could  be  used  to  Investigate  changes  In 
the  local  potentials  following  laser  Irradiation. 
Determination  of  the  effects  of  laser  irradiation 
on  choroidal  and  retinal  vasculature.  This  study 
would  be  of  clinical  Importance  In  the  use  of 
lasers  for  photocoagulatlon.  Histological 
studies  In  the  eye  and  studies  on  the  mlcroclr- 
clrculatlon  In  a simpler  system  could  be  used  to 
study  the  production  of  thromboses  In  vessels 
following  laser  exposure. 
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